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ABSTRACT 



A software package was developed to perform conceptual 
design of submarines, using the Computer v i s i on CGP-200X 
Designer System, a turn-key computer aided design 
hardware and graphics software system. The philosophy 
behind the software package is to keep all major design 
decisions under the control of the design engineer, 
rather than embedding key decisions in the program 
algor i thms. 



Modules are provided for calculating weight estimates, 
principal character istics and envelope geometry, 
resistance, weight and moment balance, and the 
equilibrium polygon. The package interfaces with a 
pressure hull design module developed separately in an 
O.E. thesis by Marvin Meade. Interactive graphics are 
used where appropriate. 



The software requires a knowledgeable naval architect as 
the user, but does not require extensive knowledge of 
computers or computer aided design systems. 
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I NTRODUCT I ON 



GENERAL DISCUSSION 



Computer Aided Design (CAD) has become an important 
tool in many industries over the past ten years. The 
capabilities of CAD systems are growing at a rate which 
defies efforts to catalog them. In the aerospace and 
automotive industries, as well as in large archi tectural 
firms, CAD has become the standard mode of design.. 

In the area of naval ship design, CAD has been 
adopted more slowly, for a number of reasons. The 
comparatively low production rates bring the issue of 
cost effectiveness into greater question for ship design 
agencies. The tremendous complexity and scale of naval 
combatants makes the design of a flexible and 
comprehensive package of software a formidable task. 

Despite these problems, a number of packages for 
conceptual and preliminary design of surface combatants 
have been implemented in both government and private 
sector design organizations. The packages used by 
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government agencies have tended to be batch-oriented 
ship synthesis programs with little or no graphic 
output » and virtually no real“t irne interactive graphic 
design capability (the notable exception to this is the 
use of interactive interior ship layout software). 

This lack of interactive graphics has been 
partially hardware driven? the concern for 
standar di zat ion and the lengthy process of government 
cer t i f i cat ion and procurement cause a substantial delay 
between the initial availability and subsequent 
installation of new hardware. In addition, the CAD 
process is foreign to many of the prominent designers in 
the very tr ad i t i on-or i en t ed wor 1 d of naval ar ch i t ectur e « 
Many design procedures in naval architecture are 
somewhat subjective and difficult to quantify, 

A further impediment to the use of CAD has been the 
tendency of many early software packages to usurp 
designer perogat ive by having key decisions in the 
design process " hard~wired" into the program algorithms. 
Such software does not utilize the expertise of the 
designer, and frequently precludes trading off design 
parameters to achieve optimal designs. 

In the area of submarine design, the above problems 
are further accentuated. The volume of new submarine 
designs is even lower than that of surface combatants. 
Accuracy becomes critical? where feet and tons may be 
acceptable units for tolerance in large surface ship 
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designs, inches and pounds are frequently the desired 
accuracy level in submarine design. This is due to the 
hydrostatic requirement to balance both surfaced and 
submerged, as well as to the premium attached to 
internal volume and deck surface area. 

Submarine designers are fewer in number than their 
surface counterparts, and they often have many design 
relationships tied to subjective criteria (i.e„, 
"designer’s eye") developed over years of experience. A 
program which does not allow the submarine designer to 
retain this subjective design flexibility is not likely 
to be accepted by the design community. 

SUBMARINE CAD: A HISTORICAL PERSPECTIVE 



Although Computer Aided Design has not been a 
primary tool in submarine design, numerous programs have 
been developed to accomplish analytical chores that are 
subsets of the overall design process. These include 
propel lor design, resistance estimates, hy drost at i c 
calculations, and structural design. These programs have 
been primarily analytical in nature, providing no 
interactive graphic capability, and are usually written 
for main-frame computer systems. 

Several packages have been written to achieve a 
complete iteration of the design process. These include 



the CODESUB program C 1 3. It provides numerical, 



rather 



-s- 



than graphic, output „ CODESUB was developed by the 
Center for Naval Analyses in order to aid in projecting 
future submarine design char ac ter i s t i cs for both our 
navy and its potential adversaries. This program was not 
intended to be a detailed tool for the naval architect, 
and, unfortunately, has many key design parameters 
imbedded in the source code, which restricts its use in 
trading off alternative designs. Other packages "graft" 
together sections of existing designs, matching the 
largest hull diameter, to provide a composite design. 
ASSET l 2 3, developed by Boeing Computer Services, while 
primarily a surface ship design package, may be expanded 
to include a submarine design module in the future. 
ASSET provides some low-level graphics output, but is 
not a full, real-time interactive graphics package. 

Thus, there is a clear opportunity for the 
development of a submarine design package incorporat ing 
the use of the full capabilities of currently available 
CAD systems. This thesis represents an exploration of 
such a design package. 

PHILOSOPHY FOR THE THESIS PROGRAM 



The software developed in this thesis accomplishes 
a "first pass" through the submarine design process at 
the conceptual design level. In conjunction with the 
concurrent Ocean Engineer thesis written by Marvin 



- 9 - 



Meade, the package includes calculation of weight 
estimates, envelope and pressure hull geometry, 
resistance (speed and power) , basic hydrostat i cs, weight 
and moment balance, and determination of the equilibrium 
polygon . 

The primary philosophy of this package is to 
develop it for a user who is a competent naval 
architect, but who may not necessarily be an expert in 
the use of computers. Crucial design decisions are left 
to the user, rather than being embedded in the program 
code. Opportunities are provided to override program 
algorithms where necessary to implement the desires of 
the individual designer. In addition, if the designer 
does posess some knowledge of CAD systems or computer 
programs, it is easy to exit the package at appropriate 
points and tailor the design to provide more detail or 
alternate geometry and analytical procedures. 
Interactive graphics are used to provide a clear 
visual i zat ion of the design under development. All 
program input parameters are couched in standard naval 
architecture terms to facilitate the user’s interaction 
with the package. The title chosen for the software 
package is CADSUB (Computer Aided Design of Submarines). 

Input to each program module is explained in the 
chapter documenting the particular module. Samples of 
program output are provided in Appendix I to this 



thesis 
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HARDUJAREs 

THE CQMPUTERV I S I ON DESIGNER SYSTEM 



GENERAL DISCUSSION 



The optimum choice of hardware for implementing the 
philosophy of this thesis was one of the several 
high-quality “turnkey" CAD systems, all of which permit 
the use of interactive graphic design techniques. 
Graphics packages for main-frame systems are available 
from several prominent suppliers, but are a compromise 
at best, since the processors they run on are not 
optimized for graphics-intensive applications. 

Suitable systems are marketed by several suppliers, 
including Computervision , Applicon, Bendix, 
McDonne 1 1 -Doug 1 as Automation Division, and others. For 
this thesis, Computervision Corporation, of Bedford, 
Massachusetts, offered free system time, instruction, 
and other technical aid to develop the program. 



THE COMPUTE RV I S I ON DESIGNER SYSTEM 



The system which Computervi si on provided time and 
instruction on is their basic “Designer System", built 
around the CGP-200X processor. The heart of the system 
consists of the processor, a high-speed tape drive, and 
one or two 300 megabyte hard disk units. The CPU is 
available with several increments of random access 
memory (RAM). The particular system used for this thesis 
was configured for appr ok i ma te 1 y 1.3 megabytes of core 
memory. The architecture of the CGP-200X is 16-bit. 

The Designer System utilizes two distinct operating 
systems. Basic file management and system-level 
housekeeping chores are handled at the “OS" level, which 
is similar in use to the operating systems on typical 
mi n i -computer installations. For the graphics operating 
system, Computervi sion uses the " CADDS 4-X" environment, 
which is tailored for graphics, and automat i cal ly 
interfaces with a powerful data base management system. 
Graphics commands are entered in a simple verb-noun 
syntax, such as "INSERT LINE" or "INSERT SPLINE". 

Input/output (I/O) is handled by several devices. 
For system administration and text file input, small 
alpha-numeric terminals are provided. For graphics I/O, 
high-resolution color or monochrome terminals are 



provided as part of a comprehensive workstation, 



which 



includes a thermal 



B.W\ d 
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printer for working quality text 
graph i cs printout. A], so provided at each work station are 
a digitising tablet and a display control device. This 
device controls background and graphics intensities, the 
number of lines of text displayed (4 or 24 lines), and 
dynamic control of the display; zooming, rotation, and 
scrolling. Output of finished quality is available from 
line printers, color pen plotters, and black and white 
e 1 ectrostat i c plotters. Large format output is supported 
by the plotters. 

PROGRAMMING LAMGU A G E . S 



Several languages are provided for use with the 
Designer System. The primary language is Fortran-S, a 
subset of ANSI Fortran. System subroutines are provided 
for interacting with the graphics and data base 
management operating systems. Programs entered in 
Fortran-S execute with the greatest speed, and allow 
creation of new graphics commands. Because of limited 
CPU space, the utilities required for this level of 
pr ogr amirri i ng , which are often provided on mia in-f ramie 
systems, must be separately programimied for each module. 
Such utilities include linking and loading capabi 1 i t ies. 
Because of the length of training required to do 
competent pr ogr amimi i ng at this level, this thesis was 



programimied in an alternate language. 
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Two macro languages are provided for the system. 
Both allow creation of execute files using standard 
graphics commands. VAR PROS is a non-comp i 1 ab 1 e macro 
language which supports Fortran-like I/O and 
computational statements. It is a flexible and effective 
language, but executes rather slowly. NEWVAR, the other 
macro language provided, is compilable, and therefore 
executes far more rapidly. The I/O for NEWVAR is 
cumbersome and difficult to format. All I/O must be 
handled in text string form, and converted within the 
program to numerical values. Despite this shortcoming, 
the execution speed of NEWVAR led to its adoption for 
use in this thesis. U1 1 imately , ^t he algorithms in this 
thesis should be translated to Fortran-S for maximum 
speed and I/O flexibility. 
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PRELIMINARY UJEIGHT ESTIMATES 



IMPORT AN CE OF WEIGHT ESTIMATES 



The algorithms employed in this thesis assume a 
weight-driven submarine design. This assumption means 
that the envelope size calculated for the submarine will 
contain enough volume to displace the same amount of 
water, when submerged, as the weight of the envelope 
displacement for the submarine. If the desired design is 
known or suspected to be volume-limited, an arbitrary 
excess of lead may be specified during the weight 
calculations, thus driving the program to calculate a 
larger envelope size. The designer may translate this 
additional weight into volume manually to allow for a 
volume-1 i mi ted design. 

Since the estimated weight data will drive the 
basic dimensions of the submarine, the weight estimating 
algorithms must be made as accurate as possible. The 
designer must carefully consider the input data to the 
weight module, and, where appropriate, may decide to 



override the programmed algorithms and specify the value 
for a particular weight group. The program is designed 
to automate the calculations, and not the design 



judgement, for the user. 


. A knowledgeable user may easily 


modify the algorithms 


themselves to suit particular 


needs - 





WEIGHT ESTIMATING ALGORITHMS 

The accounting system chosen for weight estimation 
is the SWES (Ship Work Breakdown Structure) convention 
used by the Naval Sea Systems Command for ship weight 
records C33. This system consists of seven major 
numbered groups for all ships. For submarine weight 
reports, several special weight categories are added. A 
brief summary of the weight groups is as follows: 



GROUP 


DESCRIPTION 


I 


Structura 1 


II 


Propuls i on 


1 1 1 


E 1 ectr i cal 


IV 


Combat Systems 


V 


Auxiliary Systems 


VI 


Outfit and Furnishings 


VII 


Armamen t 


A1 


Sum of I -VI I 


Lead 


Total lead 


A 


A1 + Lead 


Variable Load 


(fuel, stores, etc.) 
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NSC 



Normal Surface Condition 



Reserve Buoyancy Main Ballast- 



Submerged Disp 



Self ex planatory 



Free Flood 



Non-buoyant flooded volume 



Envelope Disp. 



Submerged Disp. + Free flood 



The algorithms chosen for the weight estimating 
module take standard design parameters as input, where 
possible, and provide an output in terms of a percentage 
of A1 weights, NSC, or submerged displacement, as 
appropriate. These percentages are then combined with 
the discrete weights that are calculated or input far 
selected weight groups to determine the A1 weight total, 
lead, NSC, and other weight parameters. The designer has 
the option of specifying weights for groups II, IV, VI, 
and VII, which override the programmed algorithms. The 
following section discusses the chosen algorithms for 
calculating each weight group. Throughout this thesis, 
the following conventions will be used to indicate 
math mat i cal o per at ions : 



"+" = addition 

= subtraction 
"-* n = multiplication 

" / " = division 
"am - exponentiation 



GROUP I : 



We i g h t 



group I is 



determined as a percentage of 
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NSC. In general, 

XW 1 =/ ( de p t h , material, size). 

The size dependence is accomodated by calculating 
W1 as a percentage of NSC. The algorithm i mp 1 emen t ed in 
the program is= 

XW1=C1+C2*DEPTH. 

Cl and C2 are dependent on hull material used, and 
are determined parametrically from data on past 
submarine designs l 4,5]. The DEPTH parameter is a user 
input of maximum operating depth in feet. The accepted 
design practice of designing for a collapse depth of 
1507. of operating depth is implemented through the 
choice of Cl and C2. 

GROUP I I s 

Weight group II is calculated directly, in tons, 
rather than as a percentage. The number of parameters 
involved in calculating W2 is greater than for any other 
weight group. 

W2=/[ horsepower, propulsion type, battery type & 
capacity (capacity^endurance) ] 

The horsepower variable for a nuclear plant is 
simply shaft horsepower, while, for a non-nuclear plant, 
two separate horsepowers are needed; shaft horsepower 
(usually the horsepower of the electric final drive 
motor) and charging horsepower (power of the prime mover 



used to drive the alternators or generators to provide 
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electrical power to recharge the storage batteries). The 
weight of the charging prime mover is a strong function 
of technology chosen; ie diesel, wankel, etc. 

The calculation of storage battery weight for a 
non-nuclear plant is dependent on battery type, capacity 
in k i lowat t-hours at a high discharge rate ("sprint"), 
and capacity in KW— H at a low discharge rate (endurance 
power). The options provided for battery type in this 
program are nickel-cadmium, improved lead-acid (German 
VARTA type), and "standard" lead acid (U.S. Trident 
type). Constants determined from the energy densities of 
these battery types at the two discharge rates are used 
to determine battery weight (the larger of sprint or 
endurance rate battery weights is chosen as the final 
battery weight). 

Electric motor weight is determined from a 
multiplier based on a modern, air-cooled DC motor. The 
option exists to provide other multipliers based on 
liquid cooling, super con due t i ng technology, etc. 

For the non-nuclear plant, an additional option 
exists to provide a lumped adjustment to weight group II 
to provide for unconventional propulsion technology 
(fuel cells, Stirling, low-power nuclear, etc.). This 
final adjustment allows the designer maximum flexibility 
to tailor group II without making modifications to the 
program code. 



For a "straight" nuclear plant, 
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W2= ( C 1 *SHP ) /[ 1 og ( SHP ) ] A C2 
where SHP=shaft horsepower 

Cl is from parametric data. 

For a non-nuclear plant, 

W2=WE+WM+WB+dW2 

where WE=charging prime mover weight 
WM=propulsion motor weight 
WB=bat tery weight 

dW2=ad justmen t to UI2 (user input) 

WE = C3* CHP 

where CHP=charging horsepower 

C3 is dependent on prime mover type. 
WM=C4*SHP. 

WB=MaxC <C5*KWHS), ( C6*KWHE ) 3 
where KWHS=capaci ty for "sprint" 
KWHE=capaci ty for endurance 
C5,C6 are dependent on battery type. 



GROUP Ills 

Although weight group III would intuitively appear 
to be a direct function of installed generating 
capacity, a study of historical design data C 5 3 reveals 
that the group III weights are, in fact, very closely 
approximated by a straight percentage of A1 weights. 
This anomaly is best explained by the large amount of 
this weight group attributable to the electrical 



distribution system spread throughout the submarine. 
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This pou jer distribution system weight is a direct 
function of the size of the submarine? since the amounts 
of cable? connectors? and distribution panels vary with 
size of the vessel. Weight group III? consequently? is 
calculated as a percentage of A1 weights: 

■/.W3=C1 . 



GROUP IV: 

Deterniination of group IV weights is very difficult 
to implement in an "automatic" mode. The variety of 
combat system components and great variation in 
equipment weights makes the selection of an algorithm 
extremely difficult? as most schemes provide accurate 
results for a limited number of combat suites? at best. 
It is therefore highly recommended that the designer 
specify a group IV weight based on his off-line 
determination of the desired combat systems and 
associated weights. To permit a rough approximation is 
acceptable for an early design iteration? the program 
contains an algorithm based on historical data as a. 
percentage of Ai weights for either fast attack or 
ballistic missile submarines: 

“/.W4 = Ci i = i ? 4 

where Ci is based on submarine type (SB? SSB? 



SSN ? or 55BN). 
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GROUP V: 

Weight group V is accurately determined as a 
function of submarine size? since the auxiliary systems 
required for a vessel are directly proportional to 
d i sp 1 a cemen t and internal volume- Group V is thus 
calculated as a percentage of A1 weights: 

7.W5-Ci i-l,2 

where Ci is based on submarine type (attack or 
ballistic missile). 

GROUP VI: 

Group VI weights are a function of crew size, and, 
to some extent, typical mission duration- The multiplier 
for group VI weight determination is based on 
habitability standards commensurate with the last two 
classes of attack and ballistic missile submarines, 
respectively- The weights are directly calculated in 
tons : 

W6=C1*NP 

where NP-number of personnel in the crew. 

GROUP VI I : 

The estimation of armament weights is nearly as 
difficult as the task of determining group IV weights, 
since the possible combinations of existing and future 
weapons and associated launchers are infinite- Again, it 
is highly recommended that the designer specify a value 



for weight group VII based on his off-line determination 
of the weapons payload. As was provided in the case of 
combat systems? a rough estimate of group VII is 
available from a percentage of A1 weights. The 
percentage constant was determined from historical data, 
and separate values are provided for attack and 
ballistic missile submarines. 

7.W7 = C 1 . 



OTHER WEIGHTS: 

The amounts of lead, variable load, reserve 
buoyancy, and free flood are provided as fractions by 
the designer. The following ranges are appropriate: 

Lead ( . 07-. 1 1 >*A1 

Variable Load ( . 04- . 07 ) *NSC (nuclear) 

( . 0S-. 19 ) *NSC ( non-nucl ear ) 



Reserve Buoyancy ( . 1 2- . 1 5 ) *NSC 

Free Flood (. 0b-. 1 )*< Submerged Displ. ) 

The operator, however, may specify his own values 
as desired. 



CALCULATION OF WEIGHT SUMMARY 



When either percentages or actual values have been 
calculated or specified for each of the weight groups, 
the calculations for the weight summary are processed. 



At a. minimum, actual weight values in tons have been 



determined for weight groups II and VI. 



Addi t 1 ona 1 1 y , 



the designer may have specified weights for weight 
groups IV and VII. As an example of the algorithms used 
for final weight calculations, assume that only W2 and 
UI6 have been determined in tons. The other permutations 
of final weight calculations are similar, but allow for 
dicrete weights, vice percentages, for either W4, W7 , or 
both. For the case where only W2 and W6 are known, the 
sequence is as follows: 

(All ' , 7. M figures used are converted to fractions; 
7.X=7.X / 1 00 ) 

K 1 = ( 7. W 1 +7.W 1 *7.LD ) / ( 1 -7. VL ) +7.W3+7.W5 
where 7L_D= lead fraction 

7.VL=var iable load fraction. 

A 1 = ( W2+W6 ) / ( 1 -K 1 -7.W4-7.W7 ) 

W3=7.W3*A1 

W4=-/.W4*A1 

W5=7.w5*A1 

W7=7.W7*A1 

Lead=7.LD*Al 

A=A1+Lead 

NSOA/< 1-7.VL) 

Var iable Load=7VL*NSC 
Main Ea 1 last=XRBY*NSC 

where 7.RBY=f ract ion of reserve buoyancy 
Submerged D i sp . =NSC+Ma i n Ea 1 1 as t 
Free FI ood=7.FF*Submer ged Disp. 



where "/.FF — free flood fraction 



Envelope Disp. -Submerged Disp. +Free Flood. 

The calculated weight summary is displayed on the 



screen with 


a summary 


of 


user inputs, 


at which time 


the 


designer may 


accept 


or 


reject the 


results. If 


the 


results are 


rejected, 


the 


program loops back to 


the 



input section for the weight module and another 
iteration is started- If the results are accepted, they 
are written into an output file for record purposes, and 
various input and calculated parameters are written into 
several intermediate files for use in later program 
modules as the design is completed. 

Because of the high degree of importance attached 
to the calculated weight data, future refinement of 
program algorithms should include a priority effort to 
accomplish two major tasks • 

(1) Better accuracy of calculated weights. 

(2) Further subdivision of weights (at least 



partially to the three-digit SWBS level). 



PRINCIPAL CHARACTERISTICS 



AND ENVELOPE GEOMETRY 

PRINCIPAL CHARACTERISTICS 

The results from the weight module provide the 
first important input to the deterrriinat ion of the 
remainder of principal character istics, which are 
geometric in nature. Since this thesis implements a 
weight-driven submarine design process, the next design 
task is to create an envelope of the appropriate volume 
to contain the envelope displacement: 

Envelope Volume=35*Envelope Displacement. 

The product envelope defined by this module of the 
program will exhibit the following geometric principal 
char act er i st i cs ( abbr evi at ions and symbols in 

parentheses are program variable names) : 

Length Overall (LOA) 

Length of Entrance (Lf) 

Length of Run (La) 

Mid Body Length (Lm) 



Diameter (D) 



Length/Diameter Ratio (L/D) 

Prismatic Coefficient (Cp> 

Bare Hull Surface Area < SF ) . 

In addition ? once envelope geometry is defined? the 
foil ow i ng calculated hydrostatic characteristics may be 
faun d : 

Submerged Longitudinal Center of Buoyancy 
( LCBSUB ) 

Draft a) Normal Surface Displacement (TNSC) 
Longitudinal Center of Buiyancy a NSC (LCBNSC). 

These parameters, then, when coupled with the input 
parameters delineated for the weight estimation module, 
constitute a set of principal characteristics for the 
des ign. 

DEFINING THE ENVELOPE 

A number of algorithms have been suggested to 
calculate the geometry of the envelope for submarines. 
These range from a " cut -an d-pas te 11 approach using past 
designs 121 to very sophisticated polynomial 
determination of offsets C6D. In between these extremes 
are a number of methods utilizing simple thumbrules and 
parametric relationships among Cp, D, L/D, and 
displacement [4,7]. 

Both the cut-and-paste and thumbrule approaches are 
somewhat imprecise in allowing the designer to fully 



control the geometry of the resulting hull form- The 
polynomial methods, while very precisely controlled, 
require the designer to spec ify design par a me t er s which 
are not couched in, or clearly related to, traditional 
naval architecture vocabulary- The most desirable 
algorithm, then, will allow strict control of the 
resulting geometry while accepting input parameters 
which are common naval archi tectural terms. 

CHOSEN ENVELOPE GEOMETRY ALGORITHM 



In choosing the algorithm implemented in this 
thesis, two restrictions were accepted for the sake of 
simplicity. The envelope shape is confined to circular 
cross-sections, and an analytically smooth envelope 
outline is calculated. Figure 4.1 demonstrates the 
nomenclature and subdivision of the envelope calculated 
with the chosen algorithm. 

The designer is prompted for the following 

geometric input parameters: 

Length of Entrance 
Length of Run 

Forebody Prismatic Coefficient (CpF) 

Aft Body Prismatic Coefficient (CpA) 

Maximum Diameter (D). 



The algorithm then defines an envelope consisting 
of a hemi—el 1 i pt i cal forebody, a cylindrical midbody, 
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and a parabolic afterbody. The chosen algorithm is a 
modification of one suggested by Jackson [4]. 




Figure 4.1: Nomenclature for determination of envelope 



geometry . 



The coefficients Nf and Na are essentially fullness 
coefficients for the fore and aft bodies, respectively. 
In order to facilitate input in traditional naval 
architectural terms, these fullness coefficients are 
calculated from the values specified by the designer for 
prismatic coefficients of the fore and aft body sections 
of the envelope: 

Nf= 71 . 0477*CpF A 5 +65 . 81 07*CpF A 4 -483 . 3037*CpF A 3 
+587. 4137*CpF A 2 -28 1 . 7224*CpF +49.5876 

Where Nf is the forward fullness coefficient. 

Na= 379. 6546*CpA A 5 -938 . 4708*CpA A 4 
+944. 8853*CpA A 3 -47 1 . 0872*CpA A 2 +1 1 9. 1 465*CpA -11.3454 
Where Na is the aft fullness coefficient. 

These polynomial relationships are valid for values 
of CpF from .5383-. 8720, and for CpA from .3333-.7111. 
These ranges subs t an t i a 1 1 y exceed those for actual 
submarine designs. The relationships were derived by 
using known ( han d-ver i f i ed ) equivalencies of prismatic 
and fullness coefficients to set up matrix equations 
with the polynomial coefficients as the unknown column 
vector. Gaussian elimination was then used to solve for 
the unkown coefficients. Additional hand-calculated 
equivalencies were then utilized to verify that, within 
the specified ranges, the polynomial relationships yield 
results correct to the fourth decimal place. 

The M X " coordinate for the overall envelope runs 



from 0 at the forward perpendicular to LOA at the aft 
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per pend i cu 1 ar . In the formula for the forebody 

hemi -el 1 i pse, a local X coordinate (Xf) runs from 0 at 
the end of the entrance length (Lf) to Lf at the forward 
per pend i cu 1 ar : 

Xf = Lf ~X 

In the aft body parabolic section, a local X 
coordinate is again used: 

Xa= X- ( Lf +Lm ) 

Note that Lm, the mid body length, is a calculated 
parameter, and not an input. The calculation of Lm will 
be described later in this chapter. 

The hull radius, r(X), for the elliptical forebody 
is calculated as follows: 

r ( X ) = R*[ < 1— ( Xf /Lf ) A Nf ) A ( 1/Nf ) ] 
where R= D/2' 

The hull radius within the mid body is simply: 
r < X ) = R 

The hull radius within the aft body parabolic 
section is: 

r ( X ) = R*C l-( Xa/La) A Na 1 

The offset r(X) is calculated for each of 21 
stations evenly spaced along the hull length. Obviously, 
r(X) at both the forward and aft perpendiculars is zero. 
When the offsets for the 21 stations have been 
determined, the sectional area of each station is 
ca 1 cu 1 ated : 

BA ( X ) = Pi*r ( X ) A 2 
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The station sectional areas are then integrated 
along the length to determine the volume enclosed by the 
en ve lope: 

rlOA 

VOL= / SA ( X ) dX 



This and all other integrations in this thesis 
done numerically, using Simpson’s Rule- 





On the first pass through this 


process, 


the 


body 


length is zero, 


and LOA= 


Lf+La. The 


resul 


volume, then, is the 


total volume 


for 


the fore 


an d 


body 


sections alone- 


To determine 


the 


pr oper 


mid 



length, the volume determined above is subtracted 
the required envelope volume, and Lm is calculated 
the length of a cylinder of radius R containing 
rema i n i ng vo 1 ume : 

DV= Envelope Volume-VOL 



are 

mid 

ting 

aft 

body 

from 

as 

the 



Lm= DV/(Pi*R A 2> 



If DV is negative, an e 
the designer is prompted for 
After Lm is calculated, 
determined : 

LOA= Lf+Lm+La 
The sectional areas for 
calculated at the new 21 
circumferences and bare 
determined : 

C i r c ( X ) = 2*Pi*r(x> 



ror message is printed, and 
new geometric input values, 
the final value of LOA is 

the resulting hull form are 
stations, after which 
hull surface area are 



r LOA 

SF= / Circ(X) dX 

J Q 

The overall prismatic coefficient is calculated by 
comparing the envelope volume with the volume of a 
cylinder (Vcyl) of length LOAs 
Vcyl= Pi*R A 2*L0A 
Cp= Envelope Volume/Vcy 1 

These calculated parameters? when taken together 
with the input values, constitute the geometric 
definition of the envelope. 

CALCULATION OF EA5IC HYDROSTATICS 

In order to later balance the submarine, it is 
important to accurately determine the longitudinal 
centers of buoyancy in the submerged (LCESUE) and normal 
surface (LCENSC) conditions. Calculation of LCENSC 
requires the determination of the draft at NSC. LCB 7 s in 
both conditions are calculated by integrating the first 
longitudinal moment of the sectional areas. 

Using the sectional areas for the fully-submerged 
hull as calculated above, LCESUE is determined: 

yrLOA 

LCBSUB= / X*SA(X) dX /VOLENV 

There are several possible approaches to 

determining the draft at NSC. When the calculations are 



done manually, it is customary to choose at least three 
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points, commonly at drafts equal to D, D/2, and D+R/2, 
and determine the volumes and resulting displacements 
for each of these drafts- A curve of displacement versus 
draft is then constructed, after which the draft 
corresponding to the NSC displacement is read from the 
curve- While this method is expeditious and acceptable 
for hand calculations, relatively large errors may 
result when the designer is not experienced enough to 
draw the proper curve through the three points- Use of 
additional points improves the accuracy, but the manual 
calculations are tedious and time-consuming. 

With the speed and accuracy of the computer 
available, a much more precise determination of draft at 
NSC (TNSC) is possible- A potentially important factor 
that is often overlooked in manual methods is the volume 
of the free flood, which is contained in the envelope, 
and may be a signifigant portion (2-107.) of the envelope 
volume. Study of existing designs indicates that 
approximately 807- of the total free flood volume is 
submerged at TNSC. Consequently, the "target" 

displacement in determining TNSC is not simply NSC, but 
rather includes this portion of the free flood: 

NSCVOL= 35*(NSC+. 8*Free Flood) 

(for 807. free flood). Where 

NSC VOL is the desired volume at TNSC. 

The method chosen to determine NSC draft is a 



bracket and halve 



iteration on the variable TNSC, 
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beginning at TNSC= 3*R/2. The allowable error in NSC 
volume is specified as 0.5 7 . in the current program. This 
tolerance may be decreased? with additional iterations 
required for convergence. The value of 0.5/1 was chosen 
as sufficient for the conceptual design phase while 
minimising the number of iterations required. The 
algorithm employed is as follows: 

DV= . 005*NSCVOL 
NPLUS= NSCVOL+DV 
NM I NUS= NSCVOL-DV 
DR= . 5*R 
TNSC= R+DR 

For each trial TNSC, the value of each station 
radius, r(X), is compared to DR: 

If ( DR ) >r ( X ) , then r(X) and SA ( X ) remain 
unchanged. 

If (DRXr(X), then: 

SA(X)= SA ( X ) -C r ( X ) A 2 *ARCOS ( DR/r ( X ) ) -DR 
* ( r ( X ) A 2 -DR A 2) A .5 D 

(Zero trim is assumed). 

/ LOA 

VOL= / SA ( X ) dX 
y 0 

If VOL>NM I NUS and VOL <N PLUS, then the correct TNSC 
is determined. If VOLCNMINUS, then DR is increased by 
. 5*DR . Similarly, if VOL>NPLUS, then DR is decreased by 
the same increment. In either case, the algorithm is 
repeated until the value of VOL is acceptable. When the 
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correct TNSC is found, LCBNSC ic calculated: 

rLQA 

LCENSO / X*SA<X) dX /NSCVOL 

7 0 

The calculated results are displayed on the screen 
for the designer’s approval. If the results are deemed 
unsatisfactory, the designer rri ay reject the envelope 
geometry and return to the input section for another 
envelope design iteration. 

If the results are accepted, all of the specified 
input values and calculated data are written to an 
output file for record purposes, and appropriate data 
for other program modules are written to passing files 
for later use. 

PLOTTING THE ENVELOPE 



An accurate plot of the envelope is imperative, 
since the outline displayed on the monitor will guide 
the designer in his interactive design of the pressure 
hull. The 21 stations determined thus far are 
insufficient to fit an accurate curve, since the 
arbitrary equal ly-spaced stations may not fall at the 
points of greatest hull curvature, or at the junctions 
of the mid body with the fore and aft body sections. 

To adequately portray hull curvature and 
transitions at section junctions, the program calculates 
132 unequally-spaced offsets to be used in plotting the 
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envelope. Offsets are most closely spaced near the 
forward perpendicular and at the transition points of 
the hull. These points are then used to fit a Bezier 
spline curve which forms the profile outline of the 
envelope. The plot is automat i cal ly scaled to make full 
use of available screen area, and both vertical and 
horizontal scales are drawn as a reference for the 
designer . 

The two-dimensional envelope outline is rotated to 
form a surface of revolution, and a mesh pattern is 
superimposed on the surface to aid visualization of the 
resulting three-dimensional hull form. In addition, an 
isometric view is displayed below the profile for 
perspect ive. 

The resulting display of two views of the outer 
hull constitutes the graphic output of the envelope 
geometry mo du 1 e . The des i gner is now ready to proceed t o 
the speed and power calculations. 
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SPEED AND POUER 



GENERAL DISCUSSION 



The calculation of resistance for the envelope 
geometry developed in the previous module is a 
straight f orwar d analytic process. The formulae employed 
are the familiar drag estimation relationships from 
basic hydrodynamics. A proposal for future 
implementation would be the addition of a propellor 
design module to accurately predict the propulsive 
coefficient, in conjunction with the envelope geometry. 

To facilitate endurance calculations for 
non-nuclear designs, the calculation of total electrical 
load at each speed is included. This total load figure 
must include all non-propul sion loads (hotel and combat 
system) as well as those required to propel the 
submar ine. 

With the data passed from the envelope geometry 
module, an estimate of appendage drag may be made. This, 



when combined with hull form and surface area data, 
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comp letes the input to the speed and power algorithm. 



ESTIMATION OF APPENDAGE DRAG 



Study of existing submarines reveals that appendage 
drag comprises thirty percent or more of the total drag 
on the submarine for most designs- In the later stages 
of the design process, detailed cons i der at i on should be 
given to the configuration and size of each appendage to 
improve the accuracy of resistance estimates- At the 
stage of conceptual design, however, details concerning 
appendages are often uncertain. Many designs are 
formulated without final selection of control surfaces, 
for example. The required shape, location and sizes of 
surfaces for cruciform, X, and inverted Y sterns may 
exhibit some variance, and the final configuration may 
not be chosen until much later in the design process. 

In view of this uncertainty, a reliable rough 
estimate of total appendage drag is desirable for the 
conceptual design- This thesis incorporates the method 
suggested by Bukalov C8H. Bukalov studied numerous 
actual submarine designs, and developed an algorithm 
based on this parametric study. This algorithm may be 
reduced to a single formula dependent on length overall 
(LOA) and maximum diameter (D). The appendage drag (Da) 
is then calculated as follows: 



Da= 1 . 09065E— 3*L0A*D+1 1 - 25 



-39- 



CHOSEN ALGORITHM 



The basic relationships i mp 1 emen ted in this program 
module may be found in numerous introductory texts 
[9,1011. The interactive input from the terminal consists 
of the following items: 

Propulsive Coefficient (PC) 

Non-propulsion Loads in KW ( HL ) 

The remainder of the required input is passed from 
the envelope geometry module, and consists of: 

Prismatic Coefficient (Cp) 

Length Overall (LOA) 

Maximum Diameter (D) 

Bare Hull Surface Area (SF) 

The speed and power data is calculated and stored 
in an array indexed on speed in knots (Vk) from zero to 
forty knots. After appendage drag is estimated in the 
manner previously described, the calculations proceed in 
an iterative loop: 

Re ( Vk ) = ( Vk* 1 . 689*L0A ) / 1 . 27908E-5 
Where Re is the Reynolds Number. 

Cf= .075/C log(Re)-2D A 2 

Where Cf is the frictional drag coefficient, 
calculated by the ITTC convention. 

Cr = Cf*C 1 . 5* < D/LOA ) A 1 . 5+7*<D/L0A ) A 3 3 



+. 002* ( Cp-. 6 ) 
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Wher e Cr is the residual, or form, drag 
coefficient, calculated after Jackson [43, and the 
last term is a correction for submarines with 
considerable parallel mid body. 

Ct = Cf +Cr+ . 00025 

Where Ct is the total drag coefficient for the 
bare hull, and the final term is the correction 
factor from tow tank studies. 

HP ( Vk > = ( 1 / PC ) *. 00872*Vk A 3*( Ct-*SF+Da ) 

Where HPCVk) is the drag in horsepower for the 
given speed, Vk. 

KW ( Vk ) = . 7457*HP ( Vk ) +HL 

Where KW(Vk) is the total load in kilowatts, 
including non-propulsion loads, for the speed, Vk. 
These calculated results for drag in horsepower and 
total load in kilowatts are displayed on the terminal 
screen for designer inspection. From these results, the 
designer may determine the maximum speed of the 

submarine for the installed shaft horsepower, as well as 
the submerged endurance at any speed for installed 
energy storage capacity. If desired, the module may be 
re-run for various values of propulsive coefficient and 
non-propul s ion loads, until a satisfactory result is 
obtained. 

When the designer accepts the calculated values, 
the results are written into an output file for 



hard-copy record purposes. 
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PRELIMINARY BALANCE 



GENERAL DISCUSSION 



When the designer has completed the speed and power 
module, he will proceed to the design of the pressure 
hull, which is a process covered separately in a 
concurrent Ocean Engineer thesis by Marvin Meade [111. 
With a graphic display of the envelope and pressure hull 
to serve as a visual aid, the designer will enter the 
subject module for this chapter to accomplish a 
longitudinal weight and buoyancy balance,, 

The purpose of the preliminary balance is to 
determine the locations of lead and main ballast tanks, 
in order to confirm the feasibility of the pressure hull 
design- The required lead (Pb) longitudinal center of 
gravity (LCG) is achieved by balancing the longitudinal 
moments of all weights except the main ballast tanks to 
create an attitude of zero trim in the normal surface 
condition (NSC). 



Following the determination of lead placement, 



an 
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add i t i ona 1 moment balance is computed to locate the main 
ballast tank LCG in order to achieve zero trim in the 
submerged condi t ion. 

In addition to the longitudinal balance, a vertical 
weight balance is computed to provide an estimate of 
transverse stability while submerged- This estimate is 
in the form of EG, the distance from the vertical center 
of buoyancy to the vertical center of gravity of the 
submarine- This parameter should have a value of about 
one foot for an optimum design. To simplify 
calculations, no trim is allowed in the surfaced 
condi t ion. 



CALCULATION PROCESS 

The first set of input data required for balancing 
the submarine is passed from previously completed 
modules. The data which is read in from earlier 
calculations is as follows: 

All 16 weight categories from the Weight Module 
( W1-W16) 

Longitudinal Centers of Buoyancy at NSC and 
submerged (LCBNSC and LCBSUB) 

Maximum Envelope Radius (R) 

Length Overall ( LOA ) - 

The first set of keyboard input is the LCG and VCG 
information for the seven primary weight groups and the 
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variable load: 

LCG ( 1 )-LCG<7> 

VCG ( 1 ) — V C G ( 7 ) 

LCG ( VL ) 

VCG ( VL ) 

Although detailed arrangements have not, as yet, 
been addressed in the design process, the designer must 
have given enough thought to rough arrangements to 
provide these location data for the aggregate weight 
subdivisions- If, when detailed arrangements are 
completed at some later time, weight group LCG and VCG 
locations are signif igant ly different from the input 
provided for this module, the balance process must be 
repeated. 

The program now calculates arrays of longitudinal 
and vertical moment data: 

LM< i ) = LCG ( i )*W< i ) 

VM< i )= VCG ( i )*W< :L ) 

Where LM and VM are the longitudinal and vertical 
moments of the weight groups, respectively. 



LCG ( Cond- A-l) = y LM(i) / y U < i ) 

i=i i=i 



7 7 _ 

VCG ( Cond- A- 1 ) = y VM ( i ) / y W(i) 

i=i i = i 



The designer is now given a screen display of the 



total weight of lead, and is asked to input the desired 
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amount of that total to be used as margin lead (ML). 
This lead will arbitrarily be located at L.OA/2 
long i tudinal ly and D/.L vertically. The remainder of the 
lead is considered to be trim lead ( TL ) . The required 
LCG for trim lead is then calculated: 

LCG ( TL ) = C W(NSC)*LCBNSC -LM(Cond. A-l ) -LM(VL) 

— ML*L0A/2 D/SL 

Similarly, the actual VCG(Lead)2 
VCG ( Lead ) = (TL*6 +ML*R > /W < Lead > 

Where 6 feet above the keel is arbitrarily chosen 
for VCG(TL) to make placement in the circular cross 
section of the ballast tanks feasible. 

The designer is then asked whether the required 
LCG(SL) is feasible. If the response is affirmative, the 
program continues on to main ballast tank (MET) 
calculations. If not, the designer is asked for new LCG 
and VCG data for the weight groups and variable load. He 
is prompted, in each case, with the previously entered 
LCG and VCG for each category. Of course, if reasonable 
adjustments in weight group locations cannot lead to 
feasible stability lead LCG (feasible being, generally, 
between the physical centers of the main ballast tanks), 
a re-appraisal of the pressure hull design and rough 
arrangements is necessary. 

When the lead has been sa t i sf actor i 1 y located, 
calculation of required MET LCG is accomplished 
(VCG(MET) is arbitrarily fixed at R)= 
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LCG(NSC)= CLM(Cond. A-l ) +LM(Lead) 

+LM ( VL ) :/W(NSC) 

LCG ( MET ) = [ Ul ( Submerged ) *LCBSUB 

— L CG ( NSC ) -*W ( NSC ) D/W ( MET ) 

The designer is again asked to decide whether the 
calculated required LCG is feasible, with similar 
results in program flow depending on the answer. 
Feasibility of MET location is dependent on the desired 
geometry of physical locations for METs as envisioned by 
the naval architect. The continued graphic display of 
envelope and pressure hull geometry is useful in aiding 
this decision process. 

The program now sums all vertical moments for 
weight groups, lead, variable load and main ballast to 
arrive at an overall VCG in the submerged condition 
(VCGSUE). An estimate of submerged stability may be 
calculated and displayed: 

EG= R-VCGSUB 

Note that the vertical center of buoyancy is 
assumed to be at R (D/2) for the submarine, which is 
basically a body of revolution. Any error due to 
buoyancy of appendages will be small and conservat i ve, 
since off-center buoyancy items (sail, etc. ) are 
generally above the geometric center of the envelope. 
Again, the designer is asked to accept or reject the 
calculated stability, and, again, rejection results in 
looping back to the input section of the program module. 
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When all calculated values are deemed acceptable, 
the input data and results are written to an output file 
for record purposes, and the locations of each item’s 
center of gravity are plotted on the pressure hull 
graphic display- The designer is now ready to proceed to 



the calculation of the equilibrium polygon- 
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THE EQUILIBRIUM POLYGON 



GENERAL DISCUSSION 

One of the most important sets of calculations in 
the submarine design process is the de ter m i na t i on of the 
equilibrium polygon. Any submarine experiences changes 
in weight distribution during the course of operations, 
as items in the variable load are consumed. A system of 
tanks must be sized and located in such a way that this 
change in weight distribution and amount may be 
compensated to allow the submarine to remain in proper 
balance and trim. 

A detailed list of the variable load items to be 
considered may be found later in this chapter. In 
aggregate, the items compose major groups in the 
categories of provisions, ammunition, stores, fuel, and 
other consumable fluids. For a nuclear or other 
non~f oss i 1 — f lie 1 submarine, the boat will always become 
lighter as items are consumed, and the compensation 



consists of adding appropriate amounts 



of variable 
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ballast at the proper locations- In the case of a diese?l 
or other fossil— fuel submarine, a paradox exists for the 
designer. As fuel is consumed from the fuel ballast 
tanks external to the pressure hull, these tanks are 
compensated with seawater, since they are "soft" tanks, 
and must remain full of fluid to maintain a zero 
differential pressure ac.ro ss the tank structure- The 
seawater used for compensation has a higher density than 
the consumed fuel which it displaces- Consequently, the 
submarine actually becomes heavier as fuel is consumed - 

The tool which is used to calculate the required 
locations and capacities of the compensating tanks is 
called the equilibrium polygon, which is a plot of 
weight versus moment (fore and aft) of the weight. 
Figure 7,1 is an example of the equilibrium polygon for 
a submarine. The points within the polygon represent 
extreme operating conditions for the submarine in terms 
of the weight and moment required to balance the 
conditions. The solid lines forming the "polygon” are 
plots of the weight and moment added by progressive 
filling and emptying of the major compensating tank 
groups. These groups, which may be comprised of several 
tanks each (the aggregate weight and moment for the 
group is plotted), are the forward trim, auxiliary, and 
after trim tank groups. If all of the operating points 
for the submarine are contained within the polygon, the 
submarine may be safely compensated for all expected 



WE IGHT IN TONS X IOO 
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FIG. 7.1 
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operat lonal conditions. If not? the designer must change 
the location(s), weight (s), or both, for the appropriate 
group ( s ) , in order to expand the polygon to include all 
the conditions. 

The conditions to be calculated for polygon points 
are specified in the NAVSHIPS Technical Manual (NSTM), 
Chapter 9290 [123. The extreme conditions are more 

severe than an actual submarine would be expected to 
encounter, and represent a "worst case" situation for 
weight and moment balance. These extreme conditions may 
be thought of as representing the following operational 
s l tuat ions: 

HEAVY # 1 • Following a short, fast patrol, during 
which no ordnance is expended, but all fuel oil is 
consumed. 



HEAVY #2: Same as the previous condition, except 

that only the fuel oil in the fuel ballast tanks is 
consumed. 

LIGHT #1 & LIGHT #2: Two variations following a 

short patrol with no fuel consumed, but with all 
or dnance expended. 

HEAVY FORWARD #1 s Following a patrol where only aft 
ordnance is expended, and all fuel oil is consumed 
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froin the forward fuel ballast tanks 

HEAVY FORWARD #2: Same as previous condition, 
except that fuel from normal fuel oil tanks is also 
consumed. 

HEAVY AFT: Fallowing a patrol with only forward 
ordnance expended, and with fuel oil from the aft 
fuel ballast and normal fuel oil tanks consumed. 

CONDITION Ns The "normal" full-load patrol 
condition at the beginning of a patrol. For 
submarines with FBTs, the FBTs are specified as 
fully ballasted with seawater. 

CONDITION Ns The same as condition N, except that 
submarines with FBTs carry a full fuel load in the 
FBTs (for such submar inss, this condition will be 
lighter than condition N>. 

For the above extreme conditions, water density is 
also specified within typical ocean ranges, so as to 
accentuate the worst case nature of thee operating 



condi t ions 



THE VAR IABLE LOAD 
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The items that comprise the category referred to as 
the variable load are n umer ous a n d d i ve r s e i n n a t ur 0 „ 
Classification of the items is again aided by the 
guidelines and categories specified in the N A VS HI PS 
Techn i cal Manual - 

The variable load is subdivided into two large 
sub-groups; fixed items, not expected to vary 
signif i cant ly during a patrol, and the truly variable 
items, routinely expended as the submarine operates. The 
"fixed” portion of the variable load consists of the 
following i terns : 

CREW AND EFFECTS: Actual weight of the submarine’s 
complement of personnel, clothing, and associated 
persona 1 effects. 

SLBMs OR COMPENSATING WATER: The weight, for a 
fleet ballistic missile (FBM) submarine, of the 
submarine launched ballistic missiles (SLBMs) or 
the compensating water in their absence. SLBM 
compensating water exactly matches the weight and 
moment of the missiles themselves, so that this 
category may, indeed, be considered a fixed item. 
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SANITARY TANKS: The weight of tanks and contents 
associated with waste disposal . 

OXYGEN CANDLES : The weight of these devices, used 
for the generation of oxygen under emergency 
conditions. These are generally a very small weight 
i tern. 

LUEE OIL IN SUMPS: That weight of lube oil 
contained in the oil sumps of equipment. 

FIXED CLEAN FUEL OIL: That weight of oil carried in 
non-com pensa ted tanks which are maintained 
essentially full (i-e., shield tanks on nuclear 
submar ines ) . 

As can be seen from Reference [i2D, there are other 



items specified as 


fixed load, 


but 


they 


are minor 


in 


nature, and beyond 


the level 


of 


detai 1 


required 


for 


conceptua 1 design ■ 


For example, 


the 


i tern 


"depth control 


tanks" in the NSTM 


is not included. 


This 


omission 


w i 1 1 



be explained in the discussion of fuel load subdivision 
in the following section. 

The items considered variable in nature are as 
foil ows : 



PROVISIONS AND STORES: The consumable foodstuffs 
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(provisions) and spare parts, paper supplies, etc. 
(stores) that the submarine crew will use during a 
patrol - 

REVITALIZATION OXYGEN: Oxygen used for life support 
(atmosphere control) during lengthy submerged 
periods. On newer submarines (particularly 
nuclear), a method of oxygen generation often 
obviates the need for oxygen banks. 

T0RPED0S, MISSILES, AND AMMUNITION: This category 
is sel f -explanatory . All expendable ordnance other 
than SLEMs is included. 

WRT TANKS': This item is the "water round torpedo" 
tank capacity, used to compensate for the presence 
or absence of torpedos in the tubes. 

RESERVE ELECTROLYTE: The electrolyte carried to 
replenish the submarine storage battery. 

FRESH WATER: For this thesis, this is a composite 
category, consisting of the total of potable, feed, 
and battery water. 

RESERVE LUBE OIL: Lube oil carried to replenish the 



sump lube oil during the patrol. 
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FUEL: Propulsion fuel for the submarine. See the 

discussion that follows for the sub— categor ies of 
fuel and their significance to the polygon 
computat ions. 

Fossil-fuel submarines present a variable load 
parodo:-; for the designer. Fuel which is carried external 
to the pressure hull (the majority of the fuel load for 
most "conventional" submarines) must be replaced by 
seawater as it is burned, since these external "fuel 
ballast tanks" (FETs) are not constructed to withstand a 
high pressure differential (i.e., they are "soft 

tanks"). The seawater which is used to compensate these 
tanks is heavier than the fuel it replaces. Thus, as the 
submarine consumes fuel, it gets heavier, rather than 
lighter. 

To provide a margin of safety at the base of the 
equilibrium polygon, some means must be provided to 
remove weight as fuel is burned. In earlier designs, a 
"safety tank" was provided, of sufficient capacity that 
when it was emptied, enough water weight was pumped 
overboard to compensate for the added seawater in the 
FBTs as fuel was consumed. The safety tank, then, was 
full at the beginning of the patrol, and would be empty, 
or nearly so, when all fuel had been removed from the 
FBTs. For purposes of this thesis, the safety tank 
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concept was not used, since the weight and volume of the 
tank and it's contents is a non-productive load item. 
However, a means to compensate for the FBT weight 
addition is necessary. 

The method chosen for implementation in this thesis 
is the provision for a variable fuel oil (VFO) tank. The 
VFu tank is internal to the pressure hull, and does not 
require seawater compensation as fuel is consumed. In 
actual practice, this tank may be ballasted and treated 
as a dual— purpose tank, acting as an additional trim or 
auxiliary tank. For computational simplicity within the 
program module, however, the VFO tank is treated as an 
uncompensated clean fuel oil tank. An operational mode 
of employment is assumed whereby the fuel in the VFO 
tank is consumed first during the patrol, before the 
FBTs are utilized. Thus, the size of the VFO tank must, 



at a minimum, 


be 


sufficient 


to compensate 


for 


the 


difference in 


fue 1 


and seawater 


densities for 


the 


FBT 



capacity. This dictates a minimum capacity of 237. of the 
total fuel load. The designer is asked to input the 
percentage of total fuel to be allocated to the VFO tank 
(or tank group, as arrangements may dictate). If 
sufficient internal volume is available, it may be 
advantageous to use a figure greater than 237., in order 
to provide an adequate margin at the base of the 



polygon . 
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CALCULAT I NG THE POLYGON 



The program module execution begins by reading the 
weight group values, submerged LCB, and LCG of the main 
ballast tanks from data files generated by the weight 
estimating and balance modules- The total load to 
submerge the submarine is calculated: 

Load to Submerg e= Submerged Displacement- Condition 
A Weights. 

The total variable load weight is displayed, and 
the designer is prompted, item by item, for the weight 
and LCG of each variable load item- As each weight and 
LCG are entered, the balance of the weight in the 
variable load account is displayed- Items which are 
commonly distributed in two sub-items, fore and aft, are 
input as two separate weights and LCGs. Examples of such 
items are torpedos, missiles, and fuel- The designer 
should have given some thought prior to this input 
sequence to the subdivision of the variable load- The 
envelope and pressure hull are displayed at the top of 
the screen during this stage of the module, to 
facilitate estimating the LCG for each item. 

The last item to be entered is fuel- The balance of 
the variable load account at this point is assumed to be 
all fuel (obviously, the balance should be zero for a 



non-fossi 1-fuel submarine) . 



The total fuel load is 
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displayed, and the designer is asked to enter the 
fraction to be allocated to the VFO tank. He is then 
asked to provide LCGs for the forward and aft FETs and 
the VFO tank. 



With the input of all items completed, the program 
displays a variable load summary of all item weights and 
LCGs. If the summary is sa t i sf ac tor y , the program 
proceeds with calculations? if not, the designer is 
returned to the input sequence, where the prompt for 
each item now includes a display of its current value. 
If this value is to remain unchanged, the operator 
simply hits the "RETURN" key? otherwise, a new value may 
be specified. This iterative input procedure is repeated 
until the designer is satisfied with the variable load 
summary . 

When the input of all variable load items is 
complete and accepted by the operator, the program 
calculates the water to balance, and the associated 
moment of that water, for each of the nine specified 
polygon points. Note that conditions M and N will be 
identical for nuclear submarines. Condition N is only 
defined for submarines with FETs. 

A computation of one of the equilibrium conditions 
would proceed in the following manner. Prior to 
computing the individual conditions, the arms and 
moments for each of the variable load items are 



determ i ned 
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Item Arm= Submerged LCB -Item LCG 
Item Moment= Item Weight *Item Arm 

In addition, factors of proportional i ty are 
determined for light and heavy seawater. For each 
condition, the sum of all variable load items, in the 
fractional amount specified by the NAVSHIPS Technical 
Manual , is computed. Similarly, the total moment of the 
specified items is determined. The required water to 
balance and moment to balance may now be calculated- 
Water to Balance= Load to Submerge -MET -Sum of 
Item Weights 

where MBT= Main Ballast Weight 

Moment to Balance- -MET Moment -Sum of Item Moments 
If required, the water and moment to balance are 
multiplied by the proport ional ity factor for light or 
heavy seawater. Each of the nine conditions is 
calculated in a similar fashion, resulting in arrays of 
weight values in tons and moments in foot-tons (with 
negative moment values indicating an aft moment). 

Several "housekeeping" chores are performed at this 
point in the module. The display of envelope and 
pressure hull, and associated data base attributes, are 
filed under the part name MIT.HULLOUT, for later recall 
to plot desired views of this graphic screen. The 
maximum values of water and moment to balance are 
examined to determine the required scale for the polygon 
plot. A new part is activated under the name MIT.POLYOUT 
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for graphically displaying the polygon. Depending on the 
scale required, one of several pre— const ructed "form" 
drawings is activated, providing the axes and 
appropriate scale information for the polygon under 
con si derat ion . 

The program now inserts and labels each of the 
equilibrium condition points on the graphic display. 
This gives the designer a clear display of the limits 
required for his polygon design. An important 
consi derat ion at this point is whether an adequate 
margin remains at the base of the polygon. Generally, if 
the lowest weight value is less than 15 tons from the 
abscissa, the designer should seriously consider 
re-structuring the variable load, as necessary, to 
achieve this margin. This may be accomplished by 
returning to the input section of the module. An example 
of such a change would be to place more of the fuel 
load, volume permitting, in the VFO tank, for a 
f oss i 1-fue 1 submar i ne . 

The operator is now asked to provide the LCGs of 
the Forward Trim, Auxiliary, and After Trim tank groups. 
He is then prompted for the capacities, in tons, for 
each of the tank groups. The program determines the arms 
and moments of each group in the same fashion previously 
delineated for the variable load items. 

Lines are now inserted on the display to represent 



the successive filling of Forward Trim, 



Aux i 1 i ary , 



and 
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After Trim tanks, followed by the successive draining of 
the tanks in the same order. These lines form the 
polygon, and all equilibrium points must fit within its 
boundaries if the tank locations and capacities are 
adequate. Obviously, a polygon whose boundaries 
substantially exceed the extent of the equilibrium 
points, while safe, represents a waste of internal 
volume, and should not be accepted. 

The designer now has the option to accept or reject 
the polygon, as plotted. If it is rejected, the polygon 
boundaries are erased, and the designer may specify new 
tank group capacities, LCGs, or both, after which a new 
polygon is plotted. This process may be repeated as many 
times as are required to achieve an acceptable polygon. 

When the polygon results are sat isf actory , the 
polygon is plotted on the e 1 ectros t a t i c plotter, the 
part MIT. POLVOUT is deleted, and the variable load 
summary, tank locations, and tank capacities are written 
to an output file for a permanent record of the final 
input data set. One iteration of the submarine 
conceptual design process, as supported by this thesis, 



is complete. 



SUMMARY AND 



RECOMMENDATIONS FOR FURTHER RESEARCH 



SUMMARY 

This thesis constitutes a first exploration of the 
use of interactive graphics for conceptual design of 
submarines. The resulting product is useful, flexible, 
and accurate for conceptual design activity. It is by no 
means a finished product in the commercial sense, but 
provides a tool to build on and to use confidently for 
such applications as student design projects. 

The power of real-time graphic display of the 
design during its creation is a major factor in the 
worth and potential of the program. Many hours of 
laborious manual drafting are avoided by the 
high-resolution graphic output. In addition, a user with 
even a rudimentary knowledge of the Computervision 
system can go far beyond the capabilities of the basic 
software package, thus greatly enhancing its value as a 



design tool 
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The program is highly interactive; the user 
controls all important design decisions- The designer 
must have a solid knowledge of naval architecture, which 
is appropriate for any ship design program- Subjective 
design judgements must remain in the province of the 
user for a credible result, and that philosophy has been 
rigorously adhered to. 

RECOMMENDATIONS FOR FURTHER RESEARCH 



There are nearly unlimited opportunities for 
continued work on this project, and in this general 
area- This thesis has scratched the surface in an area 
ripe for research and development. 

The first obvious areas encompass those portions of 
the design process which were excluded from the current 
project- These include structures, propeller design and 
detailed interior arrangements. 

Submarine structural design could benefit greatly 
from treatment in the CAD environment. Visual i cat ion of 
complex structural relationships is intrinsically 
clearer with the use of real-time graphic display. In 
addition, the use of powerful finite element packages, 
with automatic mesh generation, is possible on many CAD 
systems, including Computervision. These packages 
display stress levels in color, and will magnify 
displacements on command for rapid visual analysis. 



' 
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Ar rangemen t of interior details is feasible on 
these systems, including automatic interference checking 
and layout of piping and electrical distribution 
networks. 

Propeller design may be facilitated by graphic 
display of the flow regime during analysis. Much 
analytical work remains in this area, and research into 
links with main-frame systems and large data bases is 
also needed. 

In addition, the translation of this package into 
the Fortran-S language, with attendant programming of 
linking and loading files, would greatly improve 
execution speed and I/O flexibility. The creation of 
customized graphics commands would add more power and 
enhance the ergonomics of the package. 

Further work in the subject areas of weight 
estimation, hydrodynamics, and animated display of 
control system performance is required to derive the 
maximum benefit from submarine CAD. All of these 
subjects have been separately addressed in other 
software packages, but would be valuable additions to an 
integrated design program. 

With increased interest and design volume in the 
submarine design field, research funding is both 
necessary and clearly justified by the potential 
improvements to the speed and accuracy of the design 



process. The overall area of computer aided design of 
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submarines can support major research projects 
future? benefitting both academic and commercial 
activities. 



in t h e 
design 



— 6 6 
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WEIGHT IN TONS X IOO 
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MIT. &BCD. WEIGHT 
3-28-84 13*37144 FUTIL 6.18 



I ! WEIGHT ESTIMATE FOR 13.461 DESIGN PROJECT* 
2 ! 

31TYPE* ATTACK 

4! PROPULSION* MODIFIED WANKEL ELECTRIC 

5! battery type* nickel cadmium 

6 ! SPRI NT KW-H* 3300 
7 ! ENDURANCE KW=H* 8000 
8 1BATTERY WEIGHT * 165.631 
9!ENGINE HORSEPOWER* 4500 
JO ! ENGINE & ASSOC. EOUIP. WEIGHT * 22.86 
1.1! MOTOR & ASSOC. EOUIP. WEIGHT* 3 2.25 
1 2 ! DI SCRETE ADDITION TO GROUP 2 WEIGHT* 80 
.1 3 !M AX. OPERATING DEPTH IN FEET* 700 
14! 

1 5 ! LEAD FRACTION i 0.1 

i 6 ! variable load fraction* o.oqq 

I 7! RESERVE BUOYANCY FRACTION* 0..I25 
I 8! FREE FL(X)D FRACTION* 0.08 
19 ! 

20 ! 

21 ! CATEGORY 
22! 

23! GROUP I 
24 ! GROUP 2 
25! GROUP 3 
26 ! GROU P 4 
27 ! GROUP 5 
28! GROUP 6 
29 ! GROUP 7 
301C0ND. A I 
31 ! LEAD 
32 ! COND. A 
33! VAR. LOAD 
34 ! NSC 
35 ! MBT 

36 !SUB. DISP. 

37! FREE FLD. 

38 ! EN V. DISP. 



weight 

493.318 

300.741 

17.4418 

95 

1 37.043 

117.3 

85 

I 245.84 
124.584 
I 370.43 
1 33.884 
J 504 . 31 
1 88.039 
1692.35 
135.388 
1827.74 



M IT. &BCD. GEOMOUT 
3-28-84 !3-*38*56 FUTIL 6.18 



I iGEOMETRY OUTPUT FOR 13.461 DESIGN PROJECT* 
2 ! 

3 ! 

4! LENGTH OVERALL* 139.947 
5! MID-BODY LENGTH* 74.9465 
6! OVERALL PRISMATIC* 0.798368 
7! LENGTH/DIAMETER* 5.1832 
81SUBMERGED LCB (AFT OF FP ) * 66.2489 
9! DRAFT 0 FORMAL SURF. COND.« 21.9585 
|0!LCB 0 NORMAL SURF. COND. * 66.3329 
I I ! 

I 2! RECAP OF INPUT DATA* 

. 13 ! 

I 4 ! ENVELOPE DISPLACEMENT *1827. 74 

I 5 ! ENTRANCE* 25 

16! RUN * 40 

I 7JDIAMETER* 27 

I 8! FWD PRISMATIC* 0.6 

1 9 ! AFT PRISMATIC* 0.55 



MIT. A BCD. ()FF(HJT 
3-28-84 13i38*56 FUTIL 6.18 
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1! LENGTH WRT FP 
2!0 

310.5 
4!1. 5 
5! 2.99946 
614.99830 
716.99732 
8! 13.0646 
9 120.002 
10127.0803 
1 1134.0866 
12141.084 
13148.0813 
14155-0786 
15162.0750 
J 6 1 69.9733 
l 71 76.0706 
18183.0670 
.10100.0652 
20107.0625 
211 104.06 
221 111.057 
23! 1 18.055 
241125.052 
251132.040 
261 130. Q47 



OFFSET 

0 

1.70208 

3.4.1525 

5.08010 
6. 77120 
8. 1157 
1 1 .36/4 
13. 1284 

13.5 

13.5 

13.5 

13.5 

13.5 

13.5 

13.5 

13.5 

13.5 

13.5 

13.5 
13.334 
12.444 
10.7078 
8.07667 
4.51632 
0 
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MIT.&3CD. 
4- 9-84 



POWER 
1 0 * I3»59 



FUTIL 6. .18 



SPEED AND POWER RESULTS FOR 13.461 DESIGN PROJECT* 



non-propulsion loads (Kwm 135 
PROPULSIVE COEFFICIENT « 0.77 



7! 

8! 

9! 


SPEED <KT) 


SHP 


TOTAL KW 


1 


0.615861 


135.459 


10! 


2 


4.61783 


138.444 


1 1 ! 


3 


15.0496 


146.222 


12! 


4 


34.8408 


160.981 


13! 


5 


66.8572 


184.855 


1 4 ! 


6 


1 13.92 


219.95 


15! 


7 


1 78.818 


268.345 


16! 


8 


264.312 


332.097 


1 7! 


9 


373. 441 


413.251 


18! 


10 


508.023 


513.833 


19! 


1 1 


67 1.663 


635.859 


20! 


12 


866.75 


784 .335 


21 ! 


13 


1095.96 


952.256 


22! 


14 


1361.95 


1150.61 


23! 


15 


1 667.38 


.1378.37 


24! 


16 


2014.9 


1637.51 


25! 


1 7 


2407. 13 


1930 


26! 


18 


2846.71 


2257.79 


27! 


19 


3336.24 


2622.83 


28! 


20 


3878.34 


3027.08 


29! 


21 


4475.63 


3472.47 


30! 


22 


5130.68 


3960.95 


31 ! 


23 


5846. 1 1 


4494. 44 


32! 


24 


6624.47 


5074.87 


33! 


25 


7468.37 


5704. 1 7 


34! 


26 


8380.38 


6384.25 


35! 


27 


9363.07 


711 7.04 


36! 


28 


10419 


7904.45 


37! 


29 


11550.7 


8748.37 


38! 


30 


12760.8 


9650. 74 


39! 


31 


1405 1 .8 


10613.4 


40! 


32 


15426.3 


1 1638.4 


41 ! 


33 


16886.8 


12727.5 


42! 


34 


18435.8 


1 3882.6 


43! 


35 


20075.9 


15105.6 


44! 


36 


24809.6 


46398.4 


45! 


37 


23639.5 


1 7762.9 


46! 


38 


25568 


19201 . 1 


47! 


39 


27597.7 


207.14.6 


48! 


40 


29 731 .2 


22305.5 
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MIT. &BCD. BALOUT 
5- 6-84 10*54*3/ FUTIL 6.21 



2 1 
22 

23 

24 

25 
20 



RESULTS FROM BALANCE MODULE 



3 

A 


! GROUP 


WEIGHT 


LCG 


VCG 


5 


! 1 


660.354 


1 10 


15.5 


6 


! 2 


623.431 


120 


10 


/ 


! 3 


QQ.8 14 


1 12 


15 


a 


! 4 


105 


60 


20.5 


Q 


• 5 


243. QQ 


1 10 


Q . 5 


10 


! 6 


1Q'5.5 


10 


22.5 


1 1 


! 7 


QO 


QO 


10 


12 


! A 1 


22 1 fa . OQ 


106. 1 Q 


13.Q02 


13 


! LEAD 


221 .60Q 


76 . 88Q 


8.1414 


1 4 


1 A 


243Q.Q 


103.53 


13. 376 


15 


! V LD 


183.646 


Qfa 


16 


16 

17 


! MBT 


327. Q43 


101. 3Q 


15.5 


16 


imargin 


l LEAD (TONS ) 1 50 


0VCG=D/2 


IQ 


(stability LEAl) 


(TONS) '• 


171 . 80Q < 


20 


LIABILITY LEAD 


REQUIRED 


LCG* 67 



6 FT 



NSC (TONS)* 2623.55 ^LCB * 103.145 

SUBM (TONS)* 2Q51.4Q tfLCB * 102. Q51 
SUBMERGED STABILITY (BG): 1.72272 FT 
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MIT. OUT. &BCD. POLY 
5 . 7-B4 19t24:02 FUTIL 6.21 



1! POLYGON SUMMARY OUTPUT FILE 
2! 

3! 

4! ITEM LCG 



WEIGHT 



5 



6 ! CREW AND EFFECTS 


55 


6 


MSLBM'S 09 COMP. WATER 


0 


0 


tiiSANITARY TANKS 


85 


1.5 


D ! LUBE OIL IN SUMPS 


QO 


5 


1 0 ! FI XED CLEAN FUEL OIL 


0 


0 


1 1 ! ****TOTAL FIXED ITEMS 


**** 78.913 


11.5 


1 2! PROVISIONS AND STORES 


75 


9.5 


1 3 ! REVITALI ZAT I ON 02 


88 


1 


1 4 ! TORPEDOS IN FWD R(X)M 


45 


21 


I5!T0RPED()S IN AFT RfXJM 


0 


0 


1 6 ! TACT ICAL MISSILES FWD 


0 


0 


I7ITACTICAL MISSILES AFT 


1 16 


13 


IbiWRT TANKS 


30 


5 


ID! RESERVE ELECTROLYTE 


58 


1 


20! TOTAL FRESH WATER 


100 


1 1 


21 ! RESERVE LUBE OIL 


95 


5.5 


22! FUEL BALLAST FWD 


26.5 


28 


23! FUEL BALLAST AFT 


1 1 5 


2 1 .074 


24 ! VAR I ABLE FUEL OIL 
25! 

26! 


45 


49.074 


27 ! VAR I ABLE BALLAST TANK 


DATA : 




28! TANK GROUP 


LCG 


CAPACITY (TONS) 


2 D ! forward trim 


35 


25 


30! AUXILIARY 


75 


100 


3 1 ! AFTER TRIM 


105 


35 
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SOURCE CODE LISTINGS 



- 16 - 






VIT.^BCD. DESIGN 
5- 7—04 13:50*00 FUTIL 0.2 I 



I 'DIM W( |7) ,P( 17) .F(d) .RS(2I ) ,CIHC< 21 ) ,LX(8) ,RX(d) .XPOS(30) ,SA(2I ) <# 
2!>.LS(2I ) . H P ( 4 1 ) , KW ( 4 I ) . OU < 2 1 ) . XC 29 ) . YC 29) 

3 ! RE AD (ENTER DESIGN NAME * ) ANAME 
4!#WEIGHT CONTINUE 



5 ! N= I 
0! REPEAT 
/ ! N=N+I 

8! IF (N.E0.3.0R.N.E0.5) GOTO L(X)P 

9! PRNT DO YOU WANT TO SPECIFY <S) OR CALCULATE (C) GROUP IN] WEIGHTS? 
10! READ AH 

II! IF (&R.EQ. "S’ 1 ) F( N ) =1 

12! IF <AR.EQ."C") F(N)=0 

13! IF (F(N).EO.O) GOTO LOOP 

14! PRNT ENTER GROUP IN] WEIGHT (TONS). 

15! READ WT 

16! W(N)=WT 

I 7!#L(X)P CONTINUE 
I b! UNTIL (N. EQ. 7) 

I 9 ! READ( MAX. OPERATING DEPTH IN FEET : )ZMAX 

20 ! P ( I )=.2I55+I.606225E-4*ZMAX 

21 ! RE AD( RES EH VE BUOYANCY FRACTION: )P(I3) 

22 ! HEAD( LEAD FRACTION: )P(9> 

23 ! RE AD ( V AR I ABLE LOAD FRACTION: ) P < I I ) 

24 ! READ< FREE FLOOD FRACTION: ) P ( 15) 

25! PRNT NUCLEAR (I). DIESEL (2). OR waNKEL (3) PROPULSION? 

26! READ PM 

27! PRNT ATTACK (1) OH FBM (2>? 

2d ! READ T 

29! IF (F(2) . EQ. I ) GOTO LABI 
30 ! GOSUB GR0UP2 

3 I ! r<LAB I IF ( F ( 6 ) . EQ . I ) GOTO A I 
32!COSUB CR0UP6 

33! GOTO A I 
34 ! #GR0UP2 

35! READ (SHAFT HORSEPOWER: )SHP 

30! IF (PM. NE. I ) GOTO BATT 

37! K=.03*SHP/ I 000+ 1.13 

3d! WI=(K*IE5*SHP)/(LG(SHP )) **5 

39 ! W (2) =W 1/2240 

40 ! RTNSUB 

41 ! #B ATT CONTINUE 

42! PRNT BATrERY TYPE: NICAD (I). VARTA (2). OR TRIDENT (3)? 

43! HEAD BT 

44! READ< TOTAL ENGINE HORSEPOWER: ) ESP 

45! IF (PM.EQ.2) KE=. 0096b 

46! IF (PM.EQ.3) KE=. 0050b 

47! WE=KE*ESP 

4 o ! WW = .00645*SHP 

49! HEADtTOTAL KW-H STORAGE FOR SPRINT: ) KSP 

50! READ (TOTAL KW-H STORAGE FOR ENDURANCE RATE: ) KEND 

51 ! IF ( BT. EQ. I ) CS=35.56lO 

52! IF ( BT. EQ. 2) CS=22.353 

53! IF (BT.EQ.3) CS=I2.I925 

54! B=KSP/CS 

55! IF ( BT. EQ. I ) CE=4tt. 3 

56! IF (BT.EQ.2) CE=42.074 
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57 

56 

59 

60 
6 I 
62 
63 
04 

65 

66 
6 / 
66 

69 

70 

71 

72 

73 

/ 4 

75 
/ 6 
77 
/a 

79 
60 
d I 
62 

63 

64 

65 

66 
6/ 
66 
69 
90 
9] 

92 

93 

94 

95 

96 

97 
9d 
99 

100 
101 
102 
103 
] 0 4 

105 
1 06 
I 07 

106 
109 
I 10 
I I I 

1 1 2 
I 13 
I I 4 
I 1 5 
1 I 6 



IF (BT.b'Q.3) Cb=25 • 40 1 2 
B I =KEND/CE 
IF (Bl .GT.R) B=B 1 
WT=WE+WM+B 

PR NT ENTER DISCRETE ADJUSTMENT TO GROUP 2 (0 OR # OF TONS), 
READ DA 
W ( 2 ) =WT+DA 
RTNSUB 
#GR()UP6 

READ( NUMBER OF PERSONNEL * ) NP 
W(6>=2.3*NP 
RTNSUB 
# A 1 

IF (PM.EO. I ) P ( 3 ) -. 045 

IF ( PM. EG. 2. OR. PM.EO. 3) P(3>».0I4 

P ( 4 ) = . 06 

IF (T.EQ. I ) P ( 5 ) =. I I 
IF (T.E0.2) P ( 5 ) = . 065 
IF (T.EO. I .AND. PM.EO. 1 ) P ( 7 ) =. 03 
IF (T.EO. I .AND. PM.EO. 2) P(/)=.04 
IF (T.EO. 1 .AND. PM.EO. 3) P(7)=.05 
IF (T.EO. 2) P( /)=. 14 

FW=(P( I )+P( I ) *P(9 ) )/( 1 — P ( I I ) ) +P ( 3 ) +P ( 5 ) 

IF (F ( 4) . EQ.O . AND. F ( 7 ) . EQ.O ) GOTO CASE] 

IF (F(4) .EQ. 1 .AND.F(7) .EQ.O) G()T()CASE2 
IF (F(4) .EQ.O. AND. F(7) .EO. 1 ) GOTO CASE3 
IF ( F( 4 ) . EQ. I . AND. F ( 7 ) . EO . 1 ) GOTO CASE4 
*CASE1 

A I =(W(2)+W (6 ))/ ( 1 -FW-P ( 4 )-P( / ) ) 

W ( 4 ) =P ( 4 ) * A I 
W ( /)=P( / ) * A 1 
GOTO TOTAL 
#CASE2 

A1=(W(2)+W(4)+W(6))/(1 -FW-P ( 7 ) ) 

W( 7) =P( 7)*AI 
GOTO TOTAL 
XTASE3 

AI=(W(2)+W(6)+W (/))/(! -FW-P (4 ) ) 

W ( 4 ) =P ( 4 ) * A I 
GOTO TOTAL 
#CASE4 

Al=(W(2)+W(4)+W(6)+W( /))/(! -FW) 

*T()TAL 
W ( 3 ) =P ( 3 ) * A 1 
w ( 5 ) = P ( 5 ) * A 1 
W(d)=A ] 

W ( 9 ) =P (9 ) * A 1 
W( ] 0 ) = A 1 +W ( 9 ) 

W ( ] 2 ) =W ( 1 0 ) / ( 1 -P ( 1 1 ) ) 

W( I 1 )=P( ] | )*W( 12) 

( 1 )=P( I ) *W ( 1 2 ) 

W( I3)=P( I3)*W(I2) 

W( I 4)=W( I 2)+W( 1 3 ) 

W( I 5 ) -P ( ] 5 ) ♦W ( 1 4 ) 

W( ] 6 ) =W ( 15)+W( ]4) 

PRNT CALCULATIONS COMPLETE... 

prn r 

PRNT 

<* OUTPUT SECTION 
IF (T.EO. 1 ) *T="ATT 
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I 1 7 
lid 
I 19 
120 
121 
1 22 

123 

124 
I 25 
126 
127 
1 2d 
1 29 

130 

131 

132 

133 

134 

135 

136 

137 
1 3d 
139 
1 40 
14 I 

142 

143 
| 44 

1 4b 

146 

147 
1 4d 
1 49 
I 50 

151 

152 

153 

154 

1 55 
1 56 
15/ 

1 5d 

159 

160 
1 o I 

1 52 
1 63 

164 

165 
1 66 
167 
I 6b 
169 
1 70 
1 / I 

1 17 
1 / 3 
1 74 
I /5 

l /6 



IF (T.EQ.2) &T="BALLISTIC MISSILE" 

IF (PM.EQ.l) &PM= ,, NUCLEAH" 

IF (PM.EQ.2) &PM="DI ESEL ELECTRIC" 

IF (PM.EQ.3) &PM="WANKEL ELECTRIC" 

IF (DA.NE.O) &M="M()D IFIED " 

IF (DA.EQ.O) &M="« 

IF (BT.EQ.l) &BT= "N 1CKEL CADMIUM" 

IF (BT.EQ.2) &BT="LEAD ACID (VAWTA)" 

IF (3T.EQ.3) &BT="LEAD ACID (TRIDENT)" 

PRNT WEIGHT ESTIMATE FOR [ &N AME ) * 

PRNT 

PRNT TYPE* [ATI 

PRNT PROPULSION* [&MJC&PM] 

IF (PM.EQ. 1 ) GOTO SHAFT 
PRNT BATTERY TYPE* UBT) 

PRNT SPRINT KW-H* [ KSP ) 

PRNT ENDURANCE KW-H * [ KEND 3 

prnt battery weight* c b j 

PRNT ENGINE HORSEPOWER* [ESP) 

PRNT ENGINE WEIGHT* [WE) 

PRNT MOTOR WEIGHT* [WM] 

IF (DA.NE.O) PRNT DISCRETE ADDITION TO GROUP 2* (DA) 

#SHAFT CONTINUE 

PRNT SHAFT HORSEPOWER* [SHP) 

PRNT MAX. OPERATING DEPTH (FT)* CZMAX) 

IF ( F ( 6 ) . EO . 1 ) PRNT NUMBER OF PERSONNEL* CNP) 

PRNT LEAD FRACTION* [P(Q)J 

VARIABLE LOAD FRACTION* [ P ( II)) 

RESERVE BUOYANCY FRACTION* [ P ( I 3 ) ) 

FREE FLOOD FRACTION: [P(15)3 
HIT <RETURN > TO CONTINUE... 

ACONT 



PRNT 
PRNT 
PRNT 
PRNT 
READ 
PRNT 
PRNT 
PRNT 
1 = 1 

REPEAT 
PRNT G 
1 = 1 + 1 
UNTIL 
PRNT C 
PRNT 
PRNT 
PRNT 
PRNT 
PRNT 
PRNT 
PRNT 
PRNT 
PRNT 
PRNT 
READ 



CATEGORY 



WEIGHT 



LEAD 

COND 

VAR. 

NSC 

MBT 

SUB. 

FREE 

ENV. 



PCI) 


[W( I) ) 


EQ.d) 


i. A 1 

i 


[ W ( d) J 
t W ( Q ) 1 


i. A 


[W( 10) ] 


LOAD 


tw( 1 1 )) 
[ W ( 1 2 ) J 
[W< 13) ) 


DISP. 


[W( 14) ) 


FLD. 


[ w < 15) J 


DISP. 


tw< 1 0 ) ) 



(S) Ok RECALCULATE (R) ? 



SATISFACTORY 
&DECID 

IF ( ADECID.EQ. "R") GOTO WEIGHT 

PRNT PLEASE WAIT WHILE OUTPUT AND PASS FILES ARE WRITTEN. THE OUTPUT FOR 
PRNT THE WEIGHT CALCULATIONS WILL BE IN FILE "M IT. WEIGHT" . WHEN YOU HAVE 
PRNT COMPLETED THE ENTIRE DESIGN PROGRAM, "MIT. WEIGHT" AND THE OUTPUT 
PRNT FILES FOR SUBSEQUENT MODULES MAY BE OBTAINED IN HARD COPY BY USING 
PRNT THE COMMAND / PRINTLP"FILFNAME” / . 

OPENW 2, "MIT. WEIGHT" 

^ p — ii n 
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I 77!ASMP=SHP 
I 78 ! &ESP=ESP 
I 79!AKSP=KSP 
I t)0 ! &KEND=KEND 
I d I !&ZMAX=ZMAX 
I 82 ! &NP=NP 
I 83 ! &BB=B 
184! AWM=WM 
I 85!&HE=WE 
I86!ADA=DA 
I 87 !&LF=P (9 ) 

I 88!AVLF=P( I I ) 

I 89 ! &RBF=P ( 13) 

I 90! &FFF=P (15) 

19 | !AW|=W( | ) 

I 92!&W2=W (2) 

193! AW3=W ( 3 ) 

I 94 !&W4 =rf ( 4 ) 

I 95 ! &W5=W (5 ) 

|96!AW6=H(6) 

I 97 !&W7=W( 7) 

I 98 ! &Wd=W (8) 

1 99 ! &W9=W ( 9 ) 

200 ! &W I 0=W ( 10) 

20 1 ! &W 1 1 =rt ( | | ) 

202 ! AH 1 2=W( 12) 

203 ! &W I 3=W ( I 3) 

204 ! &W I 4=rt (|4) 

205! AW | 5=W ( 15) 

206 !&W t 6=W ( l 6) 

20 / ! &X=" WE IGHT ESTIMATE FOR "♦ANAME+ 

208 !GOSUB 0 
209 iWRITEF 2.&P 

2 I0!AX="TYPE* " +AT 
2 I 1 ! GOSUB 0 

2I2!AX=" PROPULSION* "+&M+APM 
21 3! IF (PM.EO. I) GOTO NUKE 
214! G0SU3 0 

2 i5!&x= "B attery type* "♦abt 
216! GOSUB 0 

2 I 7 ! &X="SPRI NT KW-H* "+AKSP 
2 I 8! GOSUB 0 

2 I 9 !AX = " ENDURANCE KW=H* "+AKEND 
2201GOSUB 0 

221 !AX="BATTERY WEIGHT* "+&BB 
2?2 1G0SUB 0 

22 3!3>X=" ENGINE HORSEPOWER* "+&ESP 
224 IGOSUB 0 

225 ! &X ="ENG I NE & ASSOC. EQUIP. WEIGHT* "+&WE 
226! GOSUB 0 

227 ! 4X="M()T0R A ASSOC. EQUIP. WEIGHT* " + &WM 
22a ! GOSUB 0 

229! IF (DA.EQ.J) GOTO NUKE 

230 ! 4 X="D I SCRETE ADDITION TO GROUP 2 WEIGHT* " +ADA 

231 ! GOSUB O 

232! *NUKE CONTINUE 

233 ! A X = " MAX. OPERATING DEPTH IN FEET* “♦AZMAX 
234 ! GOSUB O 
235 ! WR 1 1 'EF 2 .AP 
2 36 ! A X =" I.F Al) FRACTION* "+ALF 



237 !GOSUB () 

238!SX="VARIABLE LOAD Fri ACT ION i "♦AVLF 
239 !GOSUB 0 

240! &X=" RESERVE BUOYANCY FRACTION * " +&RBF 
24 I ! GOSUB 0 

242!&X="FREE FL(X)D FRACTION* "+&FFF 
243 ! GOSUB 0 
2 44 ! HR ITEF 2.&P 
245IWRITEF 2.AP 
240 ! &X=“CAT EGOR Y 

24 /! GOSUB O 
24d!WRITEF 2.&P 
249!&X="GR0UP I 
250! GOSUB 0 

25 I ! &X="GROUP 2 
252 ! GOSUB 0 
253!&X="GROUP 3 
254 ! GOSUB 0 
255 ! &X= JI GROUP 4 
256! GOSUB O 
257!&X="GROUP 5 
25d! GOSUB 0 
259 ! &X="GROUP 6 
260! GOSUB 0 

26 I !&X="GR()UP 7 
262! GOSUB 0 
263 ! &X="COND. A I 
264 ! GOSUB 0 
265 ! AX="LEAD 
266 !GOSUB 0 
26 / ! &X="COND. A 
268 ! GOSUB 0 
269 ! &X=" VAR . LOAD 
2 70 ! GOSUB 0 
271 ! AX ="NSC 
2 72! GOSUB O 
273 !&X="M3T 
2/4!GOSUB 0 
275 !&X="SUB. DISP. 

2/6! GOSUB 0 

2 77 ! &X = "FREE FLD. "+&HI5 

2 /d! GOSUB 0 

279!AX="ENV. DISP. "+&HI6 
2dO!GOSUB 0 
281 !GOTO DONE 
282! #0 CONTINUE 
2d3!HHlTEF 2. AX 
2d4 ! RTNSU3 
285 ! #D()NE CONTINUE 
286 ! OPENH 3 , "M IT. PASS I " 

28 / ! HR ITEF 3 . AH I 2 
288 ! WR ITEF 3.AHI5 
289! HR ITEF 3.AHI0 
290JHRI ;EF 3 .ANAME 

29 I! OPENH 4, "M IT. BALPASS" 

2921HRITEF 4. AH | 

293 ! HR ITEF 4. AH2 
294 ! WRITEF 4, &H3 
295 ! HR ITEF 4. AH4 
295 ! WR ITFF 4. AH 5 



HEIGHT" 

" + &H I 
" +&H 2 

"+&H3 

"+&H4 

"+&H5 

"+&H6 

"+&H7 

"+&Hd 

'•+&H9 

"+AH | 0 

" +AW | | 

"+&H I 2 

"+&H 1 3 



297!WRITEF 4. &W6 
298 ! WRITEF 4, AW / 
299IWRITEF 4, &W8 
300JWRITEF 4. &W9 
301 iWRITEF 4. &WI0 
3021WRITEF 4. &W| | 

303 ! WRITEF 4. AW I 2 

304! WRITEF 4, &WI3 

305JWR1TEF 4, AW|4 

306 !WR ITEF 4, AW15 

30/JWRITEF 4, SW I 6 

308 ! OP ENW 5, J 'M IT. PASSPH2 " 

309! WRITEF 5.&WI0 

3 1 0 ! WR ITEF 5.&NAME 

3 I I ! RUN NEW MIT.G3 

312! END 
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MIT.&BCD.G3 

4-27-84 I|:08i26 FUTIL 6.21 



I !D1M RS(2I ) .CIRC (2 I ) .LX ( 8 ) , RX (8 ) , XP()S( 30 ) , SA< 21 ) ,LS(2l ) <# 

2!>,HP(4I ) ,KH(4I ) ,DU(2I ) ,X( 132) ,Y( 132) . SX( 2 I ) 

3 • #GE0M 

4! OPENR I , "MIT. PASS I •' 

5 ! READF l.&W 
6 ! DNSC=&W 
7! READF l.&W 
8 ! DFF=&W 
<5! READF l.&W 
1 0 ! DENV=&W 
I I ! READF I .&NAME 
I 2! P 1=3. I4I5R265 
I3!#BRANCH CONTINUE 

I 4 ! READCENTER LENGTH of ENTRANCEl ILF 

I 5 ! READ( ENTER LENGTH OF RUNs )LA 

I6!REAU< ENTER HULL DIAMETER * ID 

I 7!READ( ENTER FOREBODY PRISMATIC COEFF. * ICPF 

I8IHEADI ENTER AFTERBODY PRISMATIC COEFF. « ICPA 

I 9! L=LF+LA 

20 ! LL=LF 

21 ! INTV=L/20 

22! R=D/2 

23!NF=7I.0477*CPF**5+65.8107*CPF**4-483.3037*CPF**3+587.4I3/*CPF**2<# 

24 ! >—28 I .7224*CPF+4Q.5876 

25!NA=379.6546*CPA**5-938.4708*CPA**4+944.8853*CPA**3-47I .0872*CPA**2<# 

26!>+l IQ. 1 465*CPA- I 1.3454 

27! LS( 11=0 

28! LS (2 1 I =L 

2Q!J=2 

30 ! REPEAT 

31 !LS(J I=( J-l )*INTV 
32! XL=LS (J) 

33 !GOSUB RAD 
34 ! RS ( J I =R ANS 
35!SA( J )-PI*RS( J I **2 
36!DU( J )=SA(J I 
37! J=J+I 

38!UNTIL( J.EQ.2I I 
3P ! GOSUB INTG 
40 ! VEN V=35*DEN V 
41 !DV=VENV-ANS 
42 ! LMB=DV/( PI *R**2 I 
43 ! LL=LF+LMB 
44 !LOA=LF+LMB+LA 
45! INTV=L0A/20 

46 ! LS (21 )=L(!A 

47 JLS< I 1=0 
48! J=2 
4Q!REPEaT 

SOILS (J)=( J-l I * I NTV 
51 ! XL=LS( J 1 
52 !GOSUB RAD 
53 ! RS ( J I =RANS 
54!CIRC(J)=2*PI*HS(JI 
55!SA(JI=PI*RS(JI**2 
56 ! DU ( .1 1 =C I RC ( J I 
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57! J=J+I 

58!UNTIL( J.EQ.21 ) 

5Q1GOSU3 INTO 
60 !SF=ANS 

6 I ! VCYL=PI*L0A»R**2 
62 !CP=VEN V/VCYL 
63! J=2 
64 ! REPEAT 

65 ! DU ( J ) =LS ( J ) *SA ( J ) 

66! J=J-H 

67!UNTIL( J.EQ.21 ) 

68 ! GOSUB INTO 
69 ! LCBSUB=ANS/VENV 
70!DR=.5*R 
71 !0N=DNSC+.8»DFF 
72! L()=DN~. 005*DN 
73!HI=DN+.005*DN 
74! PR NT DN» [ DN 1 
75 ! PRNT L()s [ LO J 
76 ! PRNT Hit CHI) 

77 ! # LCBLP CONTINUE 

78 ! PRNT ITERATING ON NSC DRAFT » T= IR+DR) 

79! J=2 
80! REPEAT 

81 ! I F (DR.GE . RS ( J > ) GOTO SAME 
82 ! AC =( AC OS ( DR/ RS ( J ) ) ) * ( P I / I 80) 

83!SX(J)=SA(J>-(RS(J> **2*AC-DR* (RS ( J ) ** 2-DR**2 ) **.5) 
84! GOTO SKIP 
85 ! #S AME CONTINUE 
86 ! SX ( J ) =SA ( J ) 

87 ! *SKI P CONTINUE 
88 ! DU ( J ) =SX ( J ) 

89 ! J =J+ I 

90 ! UNTIL ( J.EQ.21 > 

9 I !GOSUB INTG 
92!DM=ANS/35 

93! IF(DM.GE.LO.AND.DM.LE.HI ) GOTO MOM 
94! IFfDM.GT.HI ) GOTO MINUS 
95!DR=DR+. I *DR 
96 .'GOTO LCBLP 
9 7 ! #M INUS CONTINUE 
98 ! DR=DR-. I *DR 
99! GOTO LCBLP 
I 00 ! #MOM CONTINUE 
101 ! J=2 
102! REPEAT 

I 03 !DU ( J ) =LS ( J ) *SX ( J ) 

104! J=J+I 

I 05 ! UNTI L( J . EQ. 2 1 ) 

I 06! GOSUB INTG 

107! LCBNSC=ANS/ ( DM*35 ) 

I 08 ! 0X= ( I NTV-2 ) /5 
1 09 ! LX ( I )=.5 
I 1 0 ! LX ( 2 ) = 1 
I II ! LX ( 3 ) = I .5 
I I 2 ! LX ( 4 ) =2 

113! REPEAT L(X)P2t J=5. I . (J.GT.8) 

I l4!LX(J)=(J-4)*DX+2 

I 1 5 ! #L(X)P2 CONTINUE 

I I 6! REPEAT EXTRA: 1=1 .1 .(J.GT.8) 



I 
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I I 7 ! XL=LX( J ) 

I I 8 ! GOSUB RAD 

I I 0! RX( J ) =RANS 

I 20 ! #EXTRA 

121 ! GOTO GOUT 

I 22 ! #RAD CONTINUE 

I 23! IF(XL.LE.LF) GOTO FRONT 

124! IF< XL.GT.LF. AND. XL.LE. LL) RANS=R 

1 25! IF ( XL.GT. LL) GOTO AFT 

I 26! GOTO L00P4 

I 27 ! #FR0NT CONTINUE 

128! XX=LF-XL 

1 29!RANS=R* ( < I -C XX/LF) **NF) **( l/NF>> 

1 30! GOTO L(X)P4 
131 ! #AFT CONTINUE 
I 12 1 XX = XI -I I 

I 33 ! RANS=R* ( 1-( XX/LA ) **N A ) 

I 34 ! #L(X)P4 RTNSUB 

I35!#INTG <*SI MPSON'S RULE INTEGRATION 
1 36 ! ODD=0 

13 7! REPEAT ()DDL«J = 2.2. (J.GT.20) 

1 38 ! ODD=OD0+4*DU( J ) 

I 39! #ODOL CONTINUE 
I 40 ! EVEN=0 

141! REPEAT EVENL: J=3.2. (J.CT. 19) 

I 42 ! EVEN=EVEN+2*DU ( J ) 

143!*EVENL CONTINUE 

I 44! ANS=< I NTV/3 ) * ( DU ( I ) +ODD+EVEN +DU( 2 I )) 

1 45 ! RTNSUB 
I 46!#C0UT 
147! LD=LOA/D 
I 48! PRNT 
I 49 ! PRNT 

I 50! PRNT ENVELOPE DISPLACEMENT! CDENV ] 

151 ! PRNT LENGTH OVERALL: CLOAI 
1 52 ! PRNT OVERALL PRISMATIC: I CP 1 
I 53 ! PRNT DIAMETER: CD] 

I 54! PRNT LENGTH/D I AMETER : CLDJ 
I 55 ! PRNT MIDBODY LENGTH: CLMB) 

1 56 ! PRNT FWD PRISMATIC: CCPFJ 
157! PRNT AFT PRISMATIC: CCPA) 

I 58 ! PRNT ENTRANCE: CLF] 

I 59! PRNT RUN: [LAI 
I 60 ! PRNT LCB SUBMERGED: CLCBSUB] 

161! PRNT LCB 0 NSC : CLCBNSC) 

162! PRNT DRAFT 0 NSC : [R+DR] 

1 63 ! PRNT SATISFACTORY (S) OH RECALCULATE (R)7 
I 64 ! READ ADECID 

165! IF ( &DEC ID. EQ. "R") GOTO BRANCH 

1 66! PRNT PLEASE WAIT WHILE OUTPUT AND PASS FILES ARE WRITTEN. THE OUTPUT 
1 6 7 ! PRNT FILE FOR GEOMETRY SUMMARY DATA WILL BE "M IT. GEOMOUT". A FILE OF 
1 6b ! PRNT OFFSETS WILL BE IN "M IT. OFFOUT" . 

I 69 ! oPENW 3. "MIT. OFFOUT" 

1 70! &Q="LENGTH WRT FP OFFSET" 

I 7 1 ! GOSUB W 
I 72 ! LS ( I ) =0 
I 73!RS( I )=0 
1 74 !S.X=LS( I ) 

1 75 ! &C=RS ( I ) 

I 7 6 ! A S = " 
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I 77!AQ=AX+AS+AC 
I 78JWR1TEF 3 , AO 

I 79 ! REPEAT LtX)P5* J= 1 . I .(J.GT.8) 

I 80!AX=LX( J) 

181 !AC=RX(J) 

I 82 ! AQ=AX+ AS+AC 
I 83 ! GOSUB W 
I 84 ! J=J+ I 

I 85 ! #L(X)P5 CONTINUE 
I 86 ! REPEAT L(X)P6 * J =2 . 1 . ( J . GT. 2 I ) 

I 87 ! AX=LS ( J ) 

I 88 ! AC=RS( J ) 

I 8Q! AQ=AX*AS+AC 
I90!G()SUB W 
I Q I !#L(X)P6 CONTINUE 
I R2 ! GOTO FILE 
I Q3 ! #H CONTINUE 
1 94 !HR ITEF 3. AO 
195 ! RTNSUB 
I <76 ! #FI LE CONTINUE 
IQ/iOPENW 4, "MIT. PASSHP" 

1 98 ! AX=CP 

1 99 ! HR ITEF 4.&X 
200 ! AX=LOA 

20t .'WRITEF 4.&X 

202 ! AX=D 

203IWRITEF 4, AX 

204 !AX=SF 

205 !HR ITEF 4. AX 

206 ! HR ITEF 4.ANAME 

20 7 ! OPENH 5 ."MIT.GEOMOUT" 

208! AX="GE()M ETRY OUTPUT FOR "+ANAME+" **' 
209 ! WR ITEF 5 .AX 
2 10 !HR ITEF 5, AS 

2 I I iWRITEF 5 .AS 
2 I 2 ! AP=LOA 

213! &X=*'LENGTH OVERALL * " + &P 
2 I 4 iWRITEF 5,&X 
215! AP=LMB 

2 I 6 ! AX="M 1 0-BODY LENGTH: "+AP 
2 I 7! HR ITEF 5.&X 
218! AP=CP 

219! AX= J 'OVER ALL PRISMATIC* "+AP 
220!WRITEF 5. AX 
221 ! AP=LD 

222! &X=" LENGTH/O I AMETER * "+AP 
223SHRITEF 5, AX 
224 ! &P=LCBSUB 

225 ! AX="SUBMERGEO LCB (AFT OF FP)* "+AP 
226!HRITEF 5. AX 
227!TT=R+DR 
228 ! AP=TT 

229!AX="DRAFT «t NORMAL SURF. COND.* " ♦AP 
230 ! HR I TEF 5. AX 
23 1 !AP=LCBNSC 

232 ! AX="LCB a NORMAL SURF. CONO. * "+AP 

233 ! HR ITEF 5, AX 

234 ! HR I TEF 5 .AS 

235 ! A X = "RECAP OF INPUT DATA:" 

236!HRITEF 5, AX 



I 
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23 7 


WRITEF 5 • &S 




238 


&P=DEN V 




23Q 


&X= 4, ENVEL()PE 


DISPLACEMENT :"+AP 


240 


WRITEF 5,&X 




24 J 


&P=LF 




242 


&X= 4, ENTRANCE* 


"+AP 


243 


WRITEF 5 , & X 




244 


&P=LA 




245 


&X="RUN : 


J, +AP 


246 


WRITEF 5 . &X 




247 


&P=D 




248 


&X= 4 'DIAMETERS 


*'+AP 


249 


WRITEF 5.&X 




250 


&P=CPF 




251 


AX="FWD PRISMATIC: "+AP 


252 


WRITEF 5. AX 




25 3 


AP=CPA 




254 


AX="AFT PRISMATIC: "+AP 


255 


WRITEF 5, AX 




256 


OPENW 6, J, M IT. 


LCB JI 


257 


AX=LCBNSC 




258 


WRITEF 6 .AX 




259 


AX=LCBSUB 




260 


WRITEF 6. AX 




26! 


&X=R 




262 


WRITEF 6. AX 




263 


AX=LOA 




264 


WRITEF 6 , A X 




265 


OPENW 7, "MIT. 


PASSPH 1 " 


266 


AX=D 




267 


WRITEF 7. AX 




268 


AX=LOA 




269 


WRITEF 7. AX 




2 70 


AX=LF 




27! 


WRITEF /.AX 




272 


AX = LA 




273 


WRITEF 7. AX 




274 


AX=LCBSUB 




275 


WRITEF 7. AX 




2/6 


DSUB=DENV-DFF 


2 77 


&X=DSUB 




2 78 


WRITEF 7 , &X 




2 79 


& X=OENV 




280 


WRITEF /.&X 




281 


&X=NF 




282 


WRITEF 7 , & X 




283 


&X = NA 




284 


WRITEF 7 , &X 




285 


#DRAW 




286 


X ( 1 ) =0 




28 7 


Y ( ! ) =0 




288 


1=2 




289 


REPEAT 




290 


X( I ) =LX( I-i ) 




29! 


Y< I ) =RX ( I - ! ) 




292 


1=1 + 1 




293 


UNT IL( I.GT.9) 


294 


LF I NT= ( LF- ( I NTV+ 1 ))/23 


295 


XL= I NTV+ ! 




296 


I = ! 0 





2Q71REPEAT 
29d!X< I ) = XL 
2QQIG0SUB RAD 
300 1 Y ( I ) =RANS 
30 1 1 XL=XL+LFINT 
302! 1=1+1 

303 ! UNT I L ( I.GT.3 I ) 

304 ! L I NT= ( LOA-LF) / 1 0 1 

305 ! XL=LF+L I NT 

306 ! I =32 

307! REPEAT 

30d ! X ( I )=XL 

309 ! GOSUB RAD 

3 I 0 ! Y ( I ) =RANS 

3 11! XL=XL+LI NT 

312! 1=1+1 

3 1 3 !UNTIL( I.CT.I31 ) 

3 I 4 • X ( I 32 ) =L()A 
3 1 5 ! Y ( 1 32 ) =0 

3 16! OPENW d. "MIT. POINTS” 
3 I 7 ! J= I 
3 1 d ! REPEAT 
3 1 9 ! &X=X( J ) 

320 ! &Y=Y( J ) 

321 ! WRITEF 8,&X 
322!W RITEF 8,&Y 
323! J=J+1 

324!UNTIL< J.GT. 132) 

3 25 ! RUN NEW MIT. SPEED 
326 ! END 
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MIT.ABCD. SPEED 

4-11-84 I 3 « 20 * 34 FUTIL 6.18 



1 

2 

3 

4 

5 

6 

7 

8 
9 

10 
I I 
12 

13 

14 

15 

16 

17 

18 

19 

20 
21 
22 

23 

24 

25 

26 

27 

28 

29 

30 

31 

32 

33 

34 

35 

36 

37 

38 

39 

40 
4 1 

42 

43 

44 

45 

46 
4 / 

48 

49 

50 

51 

52 

53 

54 

55 

56 



DIM HP ( 40 ) . KW ( 40 ) 

READ! ENTER HOTEL AND COMBAT SYSTEM 
READ( ENTER PROPULSIVE COEFFICIENT! 
OPENR 1.-MIT. PASSHP" 
l .ACP 
1 . ALOA 
l .AD 
1 .ASF 
1 .ANAME 



LOADS 

>PC 



IN KW» ) HL 



READF 
READF 
READF 
READF 
READF 
CP= ACP 
Loa=aLoa 
D=AD 
SF=ASF 

DA= I . 09065 E-3*LOA*D+ 11.25 
J*1 

REPEAT 

RE = ( J* 1 . 689*LOA ) / 1 . 2 7908 E- 5 
CF=.075/(LG(RE)-2)**2 

CR=CF* ( 1 .5*ID/LOA) **l . 5+ 7* ( D/LOA )**3) + .002*(CP-.6) 
CT=CF+CR+. 00025 

HP ( J ) = ( l/PC)*.00872*J**3*(CT*SF+DA) 

KW(J)=.7457*HP(J)+HL 

J=J+1 

UNTIL( J.GT.40) 

PRNT 

PRNT 

PRNT 

PRNT SPEED AND POWER OUTPUT FOR I ANAME ] 

PRNT 
PRNT 
PRNT 
PRNT 
PRNT 
READ 
PRNT 

PRNT SPEED (KT> 

PRNT 
J=! 

REPEAT 
PRNT [J] 

J=J+I 

UNT I L ( J. GT. 20 ) 

PRNT HIT <RETURN> TO CONTINUE... 

READ ACONT 

PRNT SPEED (KT> SHP TOTAL KW 

PRNT 
REPEAT 
PRNT [J] 

J=J+I 

UNTIL! J.GT.40) 

PRNT PLEASE WAIT 
PRNT MODULE WILL 
OPENW 2. "MIT. POWER" 

AX="$PEED AND POWER 
WRITEF 2, AX 
AP AD I = " " 



NON PROPULSION LOADS (KW)x 
PROPULSIVE COEFFICIENT » 



[ HL J 
[PC) 



HIT <RETURN> 
ACONT 



TO SEE OUTPUT VALUES., 



SHP 



[ HP ( J ) ] 



TOTAL KW 



[ KW ( J ) ] 



[ HP ( J ) ] 



[ KW( J ) ] 



WHILE 
BE IN 



THE OUTPUT FILE IS 
FILE “MIT. POWER". 



CREATED. OUTPUT FOR THIS 



RESULTS FOR "+ANAME+ 11 * " 
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57!&PAD2=" 

58 !WR ITEF 2.&PADI 

59»&HL«HL 

60 ! APC*PC 

6J J AX*" NON- PROPULS I ON LOADS <ICW)« "+AHL 
62! WRITEF 2 ,AX 

63!AX~"P ROPULSIVE COEFFICIENT « *+APC 

64 ! WRJTEF 2, AX 
65 ! WR ITEF 2.APADI 

66»AX="SPEED (KT) SHP TOTAL KW" 

6 7 ! WR ITEF 2.APADI 

68 ! WR ITEF 2, AX 

69 ! WRITEF 2.APAD1 

70! J*1 

71 ! REPEAT 

72J&J-J 

73 ! AX*HP ( J ) 

74 ! AHP*AX( 1,8) 

75 ! AY=KW ( J ) 

76 ! AKW=AY ( 1 ,8) 

77! AOUT=* M +AJ+APAD 1 ♦ " * ♦AHP+&PAD2+AKW 

78 ! WRITEF 2.A0UT 
79! J=J+.l 

80 ! UNTI L( J . GT. 40 ) 

81 ! RUN NEW MIT.ENV 
82 ! END 



MIT. &3CD. BALANCE 
4-27-84 1 0 « 47 « 27 FUTIL 6.21 



1 ! D I M W< 16) .LC< 16) ,VC< 16) ,LM< 16) ,VM< 16) ,VY( 16) 

2 1 OPENR I .•■'MIT. BALPASS 1 ' 

3!J»I 
4 [REPEAT 
5! READF I ,&X 
6 ! W ( J ) =&X 
71 J=J+I 

8 ! UNT IL( J .GT. 16) 

0! OPENR 2. "M IT. LC3* 1 
IOIREADF 2.&X 
I I !LCBNSC=&X 
I 2! READF 2.&X 
I 3 ! LCBSUB=&X 
I 4 1 READF 2.&X 
I 5! R=&X 
t 6IREADF 2.&X 
I 7! L()A*4X 

1 8 ! PRNT THE INPUT DATA FOR THIS MODULE AND THE POLYGON MODULE WILL BE 
10! PRNT PROMPTED IN A NEW FORMAT. WHEN THE PROGRAM PROMPTS FOR AN INPUT 
20 ! PRNT VALUE, THE CURRENT VALUE WILL BE DISPLAYED TO THE RIGHT OF THE 
21 ! PRNT INPUT PROMPT MESSAGE. THIS WILL APPEAR AS « 

22! PRNT 

23 ! PRNT 'INPUT PROMPT « = CURRENT VALUE' 

24! PRNT 

25 ! PRNT THIS WILL FACILITATE CHANGES TO INPUT DATA IF MORE THAN ONE ITER- 
26! PRNT AT I ON THROUGH THE MODULE(S) IS REQUIRED. 

27!*START CONTINUE 

28! LT=0 

20 ! VT=0 

30 ! J= I 

31 ! REPEAT 

32 ! PR NT 

33! PRNT GROUP IJI LCG« 

34!READ(VALUE?)LC<J) 

35! PRNT GROUP (JI VCG« 

36!READ<VALUE?)VC(J) 

37! LM ( J ) =W ( J ) *LC ( J > 

38! VM(J )=W(J >*VC< J > 

30 ! LT=LT+LM ( J ) 

40 ! VT=VT+VM ( J ) 

41 ! J=J+ I 

42 ! UNTIL ( J . GT. 7) 

43 ! LC ( 8) =LT/W ( 8) 

44 ! VC (8 ) =VT/W (8) 

45 ! LM ( 8 ) =LT 
46! VM(8)=VT 
47 ! PRNT 

48! REAOfLCG OF VARIABLE LOAD? )LC ( I I > 

40 ! readi vcg of variable Load? > vc < 1 1 > 

50! PRNT 

5 1 ! LM ( I I > =W < || )*LC< I I ) 

52 ! VM ( I I )=W< !l )*VC< I I ) 

53! PRNT TOTAL LEAD IS [N(0)J TONS. 

54 !PRNT 

55! RE AD< AMOUNT OF LEAD FOR MARGIN IN TONS?)ML 
56!SL=W(0 >-ML 
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57 ! LSL= ( W( I2)*LCBNSC-LM(8)-LM( I I )-ML*(L0A/2 ))/SL 
58!LC(9)=(SL*LSL+ML*(L0A/2 ))/W(9) 

5P!PRNT REOUI RED STABILITY LEAD LCG IS CLSLI. 

60!READ( IS THIS FEASIBLE <Y/N>? )AH 
61 ! I F ( AR. EQ . "N" ) GOTO START 
62 ! LM ( 9 ) =W ( 9 ) *LC( 9 ) 

63!VM(9)=SL*6+ML*R 
64! VC(9)=VM(9)/W(9) 

65! LM( 10)=LM(8)+LM(9) 

66 ! VM( 1 0 ) =VM( 8 ) + VM( 9 ) 

67! LC( 10) =LM( 10)/W( 10) 

68! VC( 1 0 ) =VM ( 10)/H( 10) 

69 ! LM( 12) =LM( 10) +LM( 1 I ) 

70 ! VM ( I 2 ) =VM ( I OJ + VM ( 1 I ) 

71 ! VC ( 1 3 > =R 

72! VM ( 1 3) =R*W ( 13) 

73!LC( 13 ) = (W( I 4 ) *LCBSUB-LM ( I2))/W( 13) 

74 ! PR NT 

751PRNT REQUIRED MBT LOG IS [LC(l3)J. 

76 ! READt IS THIS FEASIBLE <Y/N>? JAR 
77! IF(AR.EQ."N») GOTO START 
78 ! LM( 1 3 ) =W ( 1 3 ) *LC( 13) 

79 ! LM ( 1 4 ) =LM ( 12)+LM( 13) 

80 ! LC ( I 4 ) =LM( I4)/W( 14) 

81 ! VM ( 1 4)=VM( 1 2)+VM (13) 

82! VC( 1 4 ) =VM( 14)/H( 14) 

83 ! BG=R-VC( 14) 

84 ! PRNT 

85 ! PRNT SUBMERGED STABILITY (BG) IS (BG) FT. 

86! READt IS THIS ACCEPTABLE <Y/N>? JAR 
87! IF(&R.EQ. "N") GOTO START 

88! PRNT THE OUTPUT FOR THIS MODULE WILL BE IN FILE *'M IT. BALOUT" . 

89 ! J= 1 

90! REPEAT 

91 !VY(J)=VC(J)-R 

92! J=J+1 

93 ! UNTIL ( J .GT. 14) 

94 ! OPENW 3, j 'MIT.BALOUT*' 

95 !&X= "RESULTS FROM BALANCE MODULE" 

961WRITEF 3,8.X 
9 7 ! & P I = " " 

98 ! &P2=" " 

9Q ! &P3=" " 

IOO!WRITEF 3. API 

101 !&X="GR()UP WEIGHT LCG VCG" 

I02IWRITEF 3. AX 
103IWRITEF 3, API 
104! J=l 
105! REPEAT 
I 06 ! AJ=J 
107 ! AX=H ( J ) 

I 08 ! AW=AX ( I .8) 

I 09 ! AX = LC ( J ) 

110! ALC=AX ( I .6) 

I 1 I ! AX=VC ( J ) 

1 I 2! AVC=AX( 1 .6) 

I I 3!A0UT=" "+AJ+" 

1 I 4 1WRITEF 3.A0UT 
I I 5 ! J =J+ I 
I 16! UNTILCJ.GT. /) 



+AW+AP2+ALC+AP3+AVC 
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I I 7 
I 18 
1 19 
120 
1 24 
1 22 

123 

124 

125 

126 
127 
.128 

129 

130 

131 

132 
1 33 
134 
J 35 

136 

137 

138 

139 
1 40 
14| 

142 

143 
1 44 
1 45 
146 
1 4 / 
148 
1 49 

150 

151 

152 

153 

I 54 

1 55 

156 

157 

158 

159 

160 
161 
162 

163 

164 

165 
1 66 
16 / 
168 

169 

170 

171 
I 12 
1 73 
1 74 
1 75 
1 /6 



J = 8 

A 1 

GOSUB 10 
J=9 

&G- " LEAD J ' 

GOSUB 10 
J= 1 0 
&G=*» A 
GOSUB 10 
J = 1 1 

&G=" V LD " 

GOSUB 10 
J=13 

&G= •* MBT 
GOSUB 10 
&NSC=W ( 12) 

&SUB =W (14) 

&LCBN=LCBNSC 

&LCBS=LCBSUB 

&BG=BG 

WRITEF 3 , &PJ 

&SL=SL 

&ML=ML 

&()UT- M MARGIN LEAD (TONS)* ML + " dVCG=D/2, LCG=L()A/2" 

WRITEF 3 ♦ &OUT 

&()UT= J, STAB I LITY LEAD (TONS)* “+&SL+” 0VCG= 6 FT 1 ' 

WRITEF 3 ,&()UT 
&LSL=LSL 

A()UT= "STABILITY LEAD REQUIRED LCG* "+&LSL 
WRITEF 3.&0UT 
W R I TEF 3 & P 1 

&OUT=-"NSC (TONS)* "♦ANSC+ J * GLCB* "+&LCBN 
WRITEF 3 ,&()UT 
WRITEF 3 *&P 1 

&()UT- J, SUBM (TONS)* "+&SUB+ J ' dLCB * , ' + &LCBS 
WRITEF 3.&0UT 
WRITEF 3.&PI 

&OUT=" SUBMERGED STABILITY (BG)* "+&BG+" FT JI 

WRITEF 3 • &OUT 

OPENW 4 , "M IT* PASS POLY" 

& X=LCBSUB 
WRITEF 4,&X 
&X=LC ( 1 1 ) 

WRITEF 4,&X 
&X=LC ( 13) 

WRITEF 4,&X 
GOTO PLOT 
#10 CONTINUE 
&X = W( J ) 

&W=&X( 1,8) 

& X = LC ( J ) 

&LC=&X ( I ,6) 

&X=VC( J) 

&VC=*X( 1 ,6) 

&OUT=&G+&W+&P2+&LC+&P3+&VC 

WRITEF 3.A0UT 

RTNSUB 

#PL()T 

SELECT MODE MODEL 
SELECT LAYER 3 
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I 77! J=l 
I 78 ! REPEAT 
I 79!&J=J 

I BO! INSERT POINTi XILC( J ) ] Y t VY ( J ) ) <CR> 

181 ! XP=LC( J ) + 1 

182! INSERT TEXT t&J) HOT I . 5« X t XP J Yt VY ( J ) I <CR> 
183 ! J=J+ I 
I84!UNTIL( J.GT. 7) 

I 85 ! I NSERT P()INT« XCLC(9) JY[VY<9) ]<CR> 

186! XP=LC(9)+I 

1 87! insert text 'LEAry hot i . 5»xtxpi yivy(Q) i<cr> 

18d! INSERT POINT* X f LC ( I 3» I Y t VY ( 13) ) <CR> 

I 89 ! XP=LC ( I 3 ) + I 

190! INSERT TEXT 'MBT' HOT I . 5 « Xt XP ] Y f VY ( I 3 ) ) <CR> 
I 9 | ! RUN NEW MIT. POLY 
I 92 ! END 



! 









m IT. ciBCD.POLY 

5- 7-84 | Q 5 1 4 so 2 FUTIL 6.21 



1 

2 

3 

4 

5 

6 
7 
d 

g 
10 
1 1 
12 
13 
1 4 
I 5 
I 6 
I / 
Id 
i g 
20 
21 
22 

23 

24 

25 

26 
27 
2d 
2Q 

30 

31 

32 

33 

34 

35 

36 

3 / 
3d 

39 

40 

4 I 

42 

43 

44 

45 

46 

47 

48 

49 

50 

51 

52 

53 

54 

55 

56 



DIM W ( I 6 ) , XV ( I 0 ) , YV ( 1 0 ) , XP ( I 0 ) , YP ( I 0) ♦ XT ( I 0 ) , YT ( 1 0 ) t M ( d ) , Y ( 8 ) 
OPENR I * J, M IT . B ALPASS J * 

J=1 

REPEAT 
RE AI)F I, AX 
W( J)=AX 
J=J+1 

UNT I L ( J . GT . 16) 

VLD=W( | | ) 

MBT=W(I3) 

DSUB=W(|4) 

OPENR 2, "MIT. PASSPOLY" 

READF 2 ,AX 
LSUB=&X 
READF 2. AX 
LVLD=AX 
READF 2 , &X 
LMBT=AX 

#5TART CONTINUE 
BAL=VLD 

PR NT TOTAL VARIABLE LOAD IS ( VLD ) TONS. 

PR NT 

RE AD( WEIGHT OF 
READ(LCG7) LCE 
BAL=BAL-CE 
GOSUB STAT 
READ ( WEIGHT OF 
READ(LCG?)LMC 

bal=bal-mc 

GOSUB STAT 
RE AD( WEIGHT OF 
RE AD( LCG? ) LST 
9aL=S AL-ST 
GOSUB STAT 
READ ( WE IGHT OF 
READ (LCG?) LLO 
BAL=BaL-L() 

GOSUB STAT 
READ (WEIGHT OF 
READ ( LCG? )LCF() 

BaL=RAL-CF() 

GOSUB STAT 
READ( WEIGHT OF 
READ ( LCG? ) LPS 
BAL=BAL-PS 
GOSUB STAT 
REA!)( WEIGHT OF 
READ (LCG?) L()2 
BAL=B AL-02 
GOSUB STAT 
RFAD(WEIGHT OF 
READ ( LCG? )LTE 
BAL=3 AL-TF 
GOSUB STAT 
READ(WEIGHT OF 
RE AD ( LCG ? ) LT A 



CREW AND EFFECTS? ) CE 



bLBM'5 OR COMPENSATING WATER? ) MC 



SANITARY TANKS AND WATER?)ST 



LUBE OIL IN SUMPS? ) LO 



FIXED CLFAN F. 0. <IE, SHIELD TANK>?)CFO 



PROVISIONS AND ST()RES?)PS 



REVITALIZATION 0XYGEN?)02 



TORPEDOS IN FORWARD RO()M?)TF 



TOEPEOOS IN AFT R (X)M ? ) T A 



- 97 - 



57 ! SAL=BaL-TA 
58 ! GOSUB STAT 

59! READtWEIGHT OF TACTICAL MISSILES FORWARD? ) MF 
00 !REaD< LCG?) LMF 

6 I ! BAL=B AL-MF 
62 JGOSUB STAT 

63! READIWEIGHT OF TACTICAL MISSILES aFT?)Ma 

64 ! READ ( LCG? ) LMA 

o5 ! BA L=3AL-MA 

661GOSUB STAT 

67!TM=TF+TA+MF+MA 

6d! LTM = (TF*LTF+.TA*LTA+MF*LMF+MA*LMA)/TM 
69 ! READ ( WE IGHT OF WATER FOR WHT TANKS7JWRT 
/0!READ(LCG?)LWRT 

7 I !BaL=BaL-WRT 
/2!G()SUB STAT 

73 ! RE AD ( WE IGHT OF RESERVE ELECTROLYTE?) EL 
74 ! READ t LCG?) LEL 
75 ! B AL=BaL-EL 
/6 SGOSUB STAT 

77! READ! WEIGHT OF FRESH WATER <P()TT ABLE .FEED . BATTERY RESERVE>?)FW 
7a ! READ! LCG?) LFW 
/9 ! BAL=3 AL-FW 
aO!GOSUB STAT 

8 I ! RE AD (WEIGHT OF RESERVE LUBE & HYDRAULIC OIL?)RL() 

82 ! READ( LCG? ) LRLO 

83! BaL=BAL-RL() 

84 ! GOS'JB STAT 

85 ! PRNT THE LAST VARIABLE LOAD ITEM TO BE CONSIDERED IS FUEL OIL. TO ENSURE 
86 ! PRNT A MARGIN OF SAFETY AT THE BASE OF THE POLYGON. A VARIABLE FUEL OIL 
B7 ! PRNT ( VFO) TANK IS AUTOMATICALLY SPECIFIED. THE REQUIRED VFO CAPACITY IS 
88! PRNT COMPUTED USING YOUR INPUT OF THE FRACTION OF TOTAL FUEL TO 
89 ! PRNT PLACE INTO THE VFO TANK. THE MINIMUM SUGGESTED FRACTION IS 
90 ! PRNT .23 (THE RATIO DIFFERENTIAL OF FUEL AND WATER DENSITIES). 

9 | ! PRNT IT IS HIGHLY RECOMMENDED THAT AN ADDITIONAL MAHGIN OF AT 
92 ! PRNT LEAST .05 BE SPECIFIED TO ALLOW A SAFETY MARGIN AT THE BASE 
93! PRNT OF THE POLYGON. 

94! PRNT 

95: re ad( fraction of fuel to be placed in vfo?> vfr 

96!TF()=BAL 

9 /! VFO=VFR*TF() 

98 !FBT=TF()-VF() 

99! FL=CE+MC+ST/3+L0+CF0 

I 00!LFL=(CE*LCE+MC*LMC+ST*LST+L()*LL()+CF()*LCF0)/FL 
101 !PRNT FBT CAPACITY (TONS): IFBTI 
102! PRNT VFO CAPACITY (TONS): IVFOJ 
I 03 ! PRNT 

1 04 ! RE AD( AMOUNT OF FBT <T()NS> IN FORWARD TANKS?) FBTF 

I05!FBTA=FBT-FBTF 

106! PRNT 

IO/!READ(FWD FBT LCG71LFBTF 
1 08! HEAD ( AFT F3T LCG?)LF3TA 
1 09 ! LFBT= ( FBTF*LFBTF+FBT A* LFBT A ) /FBT 
I 10! RE AD (VFO LCG7JLVF0 
III! PRNT 
I I 2 ! PRNT 
I I 3 ! PRNT 
I I 4 ! PRN r 

I I 6 ! PRNT VARIABLE LOAD RECAP: 

I 16! PRNT ITEM WEIGHT LCG 
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1 1 / 



116! 


PRNT 


CREW & EFFECTS 


(CEJ 


[ LCE ] 


1 IQ! 


PRNT 


SLBM'S 


[MCI 


[ LMC ] 


120! 


PRNT 


SANITARY TANKS 


(ST) 


[ LSTJ 


121 ! 


PRNT 


SUMP LUBE OIL 


(LO) 


( LLO] 


122! 


PRNT 


CLEAN F. ()• (FIXED) 


[CFO) 


[ LCFO ] 


123! 


PRNT 


***T()TAL fixed*** 


(FL) 


f LFLI 


1 241 


PRNT 


PROV. & STORES 


(PS I 


[ LPS ) 


125! 


PRNT 


REVITALIZATION 02 


(021 


C L()2 I 


1 26! 


PRNT 


TORPEDOS FWD ROOM 


(TF) 


f LTF ] 


127! 


PRNT 


TORPEDOS AFT ROOM 


[TAJ 


[ LTA ] 


1 28! 


PRNT 


TACT. MISSILES FWD 


(MFJ 


[LMFJ 


12Q! 


PRNT 


TACT. MISSILES AFT 


(MAI 


[ LM A] 


1 30! 


PRNT 


WRT TANKS 


(WRT) 


[LWRT) 


131 ! 


PRNT 


RESERVE ELECTROLYTE 


(ELI 


[LELJ 


1 32! 


PRNT 


TOTAL FRESH WATER 


[FW] 


[ LFW ) 


1 33! 


PRNT 


RESERVE LUBE OIL 


[ RLO] 


[ LRLO] 


134! 


PRNT 


FUEL BALLAST FWD 


[ FBTF) 


[ LFBTF J 


135! 


PRNT 


FUEL BALLAST AFT 


[ FBTA ] 


[ LF3TA ] 


136! 


PRNT 


VARIABLE FUEL OIL 


[ VFOJ 


[ LVFO ] 



138 
1 IQ 
1 40 
1 4 1 
I 42 
I 43 
I 44 
145 
1 46 
147 
1 46 
I 4Q 
I 50 
151 
1 52 
1 53 
154 
1 55 

156 

157 

158 
I 5Q 
160 
161 
162 
163 
1 64 
1 65 
I 60 
167 
1 66 
I 6Q 
1 70 
I l\ 

I 72 
I / 3 
1 74 
1 / 5 
I / 6 



PRNT 



I F( AR . EQ . "NO” ) GOTO START 
PRNT 

AFL=LSUB-LFL 

APS=LSUB-LPS 

A02=LSUB-L02 

ATM=LSUB-LTM 

AWRT=LSUB-LWRT 

AEL=LSUB-LEL 

AFW=LSUB-LFW 

ARLO=LSUB-LRL(> 

AFBT=LSUB-LFBT 

AVFO=LSUB-LVFO 

aTF=LSUB-LTF 

ATA=LSUB-LTA 

AMF=LSUB-LMF 

ama=lsub-lma 

AFBTF=LSUB-LFBTF 

AF BTA=LSUB-LFBT A 

MFL=FL*AFL 

MPS =PS* APS 

M02=()2*A02 

mtm=tm*atm 

MWRT=WRT*AWRT 

MEL»EL*AEL 

mfw«fw*afw 

MRL0=RL()*ARL<) 

MFBT=FBT*AFBT 

MVF()=VF()*AVFO 

MTF=TF*ATF 

MTA=TA*ATA 

mm.f=mf*amf 

WMA=Ma*AMA 
MFBTF=FBTF*AFBTF 
MFBTA=FBTA*AF3TA 
CVL= . QQ3 75 
CVH=1 • 00Q3 /5 
GOTO NSTM 
*5TaT CONTINUE 
PRNT 
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17/! PRNT BALANCE ()F VARIABLE LOAD ACCOUNTx [BAL] 

I 7B ! PRNT 
I /Q ! RTNSUB 
1 80 ! #NSTM CONTINUE 
Idl !LDSUB=W( 1 4 ) -H (10) 

1 82 ! #HEAV 1 CONTINUE 
l83!CW=43/35 

184!HI=FL+TM+.5*PS+EL+WRT+.5*RL()+.5*FW+CW*FBT 
I 85 ! WH 1 =MFL+MTM+. 5*MPS+MEL+MWRT+. 5*MRL(>+. 5*MFW+CW*MFBT 
I 86 ! VBH 1 =LDSUB-MBT-H I 
I 87 ! AMBT=L5UB-LM3T 
I 88! MUBT=MBT*AMBT 
I 80 ! MVBHI =-MMBT-MH1 
I 00! XV ( I ) =MV3H I *CVL 
1 0 I ! YV ( 1 )=VBH1 *CVL' 

102! #HEAV2 CONTINUE 

I03!H2=FL+TM+PS+EL+HRT+RLO+FH+CW*FBT 

IQ4!MH2=MFL+MTM+MPS+MEL+‘JlWRT+MRLO+MFW+CW*MFBT 

I 05 ! VBH2=LDSUB-MBT-H2 

I06!MVBH2=-MMBT-MH2 

107! XV (2) =MVBH2*CVL 

|08!YV(2)=VBH2*CVL 

1 00 ! #L ITE I 

200 ! L I =FL+. 75*PS+WRT+. 75*RL(H-. /5*FH+FBT+VF(> 

20.1 !MLI=MFL+.75*MPS+MWHT+. /5*MRL()+. /5*MF*+MFBT+MVF0 

202 ! VBL I =LDSUB-MBT-L I 

203! MVBL1 =-MMBT-MLI 

204! XV(3)=MVBLI *CVH 

205 ! YV ( 3) =VBL1 *CVH 

206!#LITE2 

20 7 ! L2=FL+WRT-t-RL()+VF()+CV<*FBT 

208 ! ML2=MFL+MWRT+MRL()+M VFO+CW*MFBT 

200!VBL2=LDSUB-MBT-L2 

210!MVBL2=-MMBT-BL2 

2 II ! XV ( 4 ) = MVBL2*CVH 
212!YV(4)=VBL2*CVH 
213! #HF 1 

214!HFI=FL+MF+TF+. 75*PB+EL+WRT+. 75*RL(>+. 75*FW+VF()*FBTA+CW*FBTF 

2I5!MHF1»MFL+MTF+MUF+. /5*MPS+MEL+MWRT+. /5*MRLc>+. /5*WFH+MVF()+MFBTA+CW*PFBTF 

216! VBHF1 =LDSUB-MBT-HF1 

21 7 ! MVBHF I =-MMBT-MHFl 

2 1 8 ! XV ( 5) =MVBHF I 

2 1 0 ! YV ( 5 ) =VBHF I 

220!#HF2 

72 1 !HF2=FL+MF+TF+. 5*Po+EL+WRT+. 25*HL()+. 5*FW+VF()+FBTA+CW*FBTF 

222!MHF2=MFL+MTF+MMF+.75*MPS+MEL+wwRT+.25*MRLO+.5*»/FW+MVFO-t-MFBTA+CW*MFBTF 

723 ! VBHF2=LD8UB-MBT-HF2 

224 ! M VBHF2=-MMBT-MHF2 

725!XV(6)=MVBHF2 

226! YV (6 ) =VBHF2 

227! #HA 

728 ! HA =FL+MA+TA+. 5*PS+*RT+ . 75 *HL<)+ . 5*FW +VF()+FBTF+CW*FBTA 

220!MHA=MFL+MMA+MTA+.5*MPS+MWRT+. 75*«RLo+. 5*MFWf VVFO+MFBTF+CW»MFBTA 

730! VBHA=LI)SUB-MBT-HA 

231 ! VVBHA=-MMBT-MHA 

232!XV( 7)=MVBHA 

233 ! YV ( / ) =VBHA 

234 ! #N 

735! N=H2 

236!MN=WH2 



23 / 
23d 

239 

240 

24 I 

242 

243 
2 44 

245 

246 

247 
24d 

249 

250 

251 

252 

253 

254 

255 

256 

257 
25d 

259 

260 
261 
262 

263 

264 

265 

266 
267 
26d 
269 
2/0 
271 
2/2 

273 

274 
2/5 

276 

277 
2 76 
2 /9 
280 
281 
282 

283 

284 

285 

286 
287 
2 8d 

289 

290 
29 l 
29 2 
29 3 
29 4 

295 

296 






FORWARD TRIM TANK GROUP7JLFT 
AFTER TRIM TANK GROUP7JLAT 
AUX TANK GR0UP7JLAUX 



V3N=LDSUB-MBT-N 
MVBN=- MMBT-MN 
XV (8 ) =MVBN 
YV ( 8 > =VBN 
#CONDM CONTINUE 
CM=H2-CW*FBT+FBT 
MCM=MH2-CW*MFBT+MFBT 
VBCM=LDSUB-MBT-CM 
MVBCM=-MMBT-MCM 
XV ( 9 ) =MVBCM 
YV ( 9 > =VBCM 

<#***★★****★★** SET UP LOCATION OF TRIM * AUX TANK'S 
PRNT 

PRNT YOUR ACTUAL DESIGN MAY HAVE AS MANY TRIM AND AUX TANKS AS YOU DEEM 
PRNT APPROPRIATE FOR YOUR PARTICULAR ARRANGEMENT SCHEME. TO SIMPLIFY THE 
PRNT CALCULATION PROCESS FOR THE EQUILIBRIUM POLYGON , YOU ARE ASKED TO 
PRNT ENTER AGGREGATE CAPACITY AND LCG FOR THE FORWARD TRIM GROUP. THE 
PRNT AFTER TRIM GROUP, AND THE AUXILIARY GROUP TANKS. 

PRNT 

GOTO SIZE 

#L(X)P 

EXECV 

READ (LCG OF 
READ( LCG OF 
READCLCG OF 

aft=lsub-lft 
aat=lsub-lat 

AAUX=LSUB-LAUX 
RTNSUB 

#S I ZE CONTINUE 
J=1 

REPEAT 

I F ( XV ( J ) . GT . 5 000 > GOT 0 LARGE 
IF (XV (J ) .LT. -5000) GOTO LARGE 
I F ( YV ( J ) .GT.500) GOTO LARGE 
J=J+ 1 

UNTIL ( J • GT • 9 ) 

GOTO SMALL 

# large 

DEL PAR MIT. HULLOUT 
EXIT PART FILE MIT. HULLOUT 
XS=1000 
Y S = 1 00 

ACT PAR M IT. POLYOUT 

ACT DRA POLYGON FORM M IT . LaRGEPOLY DRAW POLY 
SEL CPL LEFT 

DEF VIE GRAPH: X l l Y6 , XOYO , X22Y 1 / <CR> 

GOTO GRAPH 
45MALL CONTINUE 
DEL PAR MIT.HULLOUT 
EXIT PART FILE MIT.HULLOUT 
XS=500 
YS=50 
DEL PAR 
ACT PAR 
ACT DRA 
SEL CPL 
DEF VIE 



#GRAPH 

J=1 



MIT. POLYOUT 
MIT. POLYOUT 

POLYGON FORM M IT. SM A LLPOLY DRAW POLY 2 
LEFT 

GRAPH: X 11 Y6.XOYO,XP2Y 1 / <CR> 

CONTINUE 
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298 

299 

300 

301 

302 

303 

304 

305 

306 

307 

308 

309 
3 10 
31 1 
3 12 

313 

314 

315 

316 

317 

318 

319 

320 

321 
3 22 

323 

324 

325 

326 

327 

328 

329 

330 
33 1 

332 

333 

334 

335 

336 

33 7 

338 

339 

340 

341 

342 

34 3 
3 44 

345 

346 
34 7 

343 

349 

350 

351 

352 

353 

354 
3 55 
356 



REPEAT 

XP(J)=XV(J)/XS 
YP( J)=YV( J )/YS 

INSERT POINT: XI XP ( J ) ] Y I YP ( J ) 1 <CR> 

XT ( J ) =XP( J ) ♦. 2 
YT( J) =YP( J )-. 1 
J-J + 1 

UNTIL(J.GT.Q) 

INS TEX HI HG .2: XI XT ( 1 ) ) Y I YT ( 1 ) ] <CR> 

INS TEX H2 HG .2: X I XT ( 2 ) J Y I YT ( 2 ) J <CR> 

INS TEX LI HG .2: XI XT (3 ) ] Y I YT ( 3 ) ) <CR> 

INS TEX L2 HG .2: X I XT( 4 ) ] Y I YT( 4 ) ] <CR> 

INS TEX HF1 HG .2: X I XT ( 5 ) 1 Y I YT< 5 ) ] <CR> 

INS TEX HF2 HG ,2: X I XT(6 ) ] Y [ YT( 6 ) ] <CR> 

INS TEX HA HG .2: X I XT ( 7 ) 1 Y I YT ( 7 ) ] <CR> 

INS TEX N HG .2: X I XT ( 8 ) - . 5 J Y I YT ( 8 ) ] <CR> 

INS TEX M HG .2: XI XT ( 9 ) ] YI YT(9 ) 3 <CR> 

gosub Loop 

GOTO PLTNK 
#C AP continue 
EXECV 

READICAPACITY <TONS> of Forward TRIM GROUP?) CFT 
READICAPACITY <T()NS> OF AUXILIARY GR()UP?)CAUX 
READICAPACITY <T()NS> OF AFTER TRIM GROUP?)CAT 
RTNSUB 

#PLTNK CONTINUE 
GOSUB CAP 
M ( 1 ) =0 

M(2)*AFT*CFT 

M (3 )=M(2 ) + AAUX*CAUX 

M(4)=M(3)+AAT*CAT 

M (5 > =M( 4 ) -A FT* CFT 

M(6)=M(5)-AAUX*CAUX 

M ( 7 ) =0 

Y ( 1 )=0 

Y ( 2 ) =CFT 

Y ( 3 ) =CAU X+CFT 
Y(4)=Y(3)+CAT 

Y ( 5 ) =Y ( 4 )-CFT 
Y ( 6 ) =Y ( 5 ) -CAUX 

Y ( / ) =0 
J=l 

REPEAT 

M ( J ) =M ( J )/ XS 
Y( J)=Y( J)/YS 
J = J+1 

UNTIL! J.GT. 7) 

INS LIN TAG=1 : X I M ( I ) ]YIY(I ) J * X [ M ( 2 ) ] Y [ Y ( 2 ) ] « <CR> 

X I M ( 2 ) 1 Y I Y ( 2 ) ) , X I M ( 3 ) J YI Y ( 3 ) J ♦ <CR> 

X I M ( 3 ) ] Y I Y (3) I * X I M ( 4 ) JYIY 14) ] , <CR> 

INS LIN TAG=4i X [ M ( 4 ) J Y C Y ( 4 ) ] • X I M ( 5 ) ] Y I Y ( 5 ) ) . <CR> 

INS LIN TAG=5: X I M (5 ) J Y I Y ( 5 > J , X I M (6 ) ] Y I Y ( 6 ) J , <CR> 

INS LIN TAG=6 : XI M ( A ) ] YI Y ( 6 ) ] • X f M( 7 ) ] Y[ Y( 7 ) ] , <CR> 

EXECV 

READ! IS THIS POLYGON SATISFACTORY <Y/N> ? ) A R 

IF! AH.EQ. ,, Y 1 ' ) GOTO FINAL 

DEL ENT: TAG 1 .TAG 2. TAG 3. TAG 4. TAG 5. TAG 6 <CR> 

EXECV 

REA!) (CHANGE CAPACITIES ONLY <1>. OR LOG'S & CAPACITIES <2> ? 
I F ( CHC. EO . 2 ) GOTO BOTH 



INS LIN TA 0=2: 
INS LIN TAG=3: 



) CHC 



357! IF(CHC.E0. I ) GOTO PLTNK 

358 ! #B()TH COtTTlNUE 

35 < ?!G0SUB L(X)P 

300 ! GOTO PLTNK 

361 !#FINAL CONTINUE 

362 ! OPENW 3 , J, M IT. PASS. POLYDAT" 

363 ! & X=CE 

364 ! GOSUB 0 

365 ! &X=LCE 

366! GOSUB 0 

367 !&X=MC 

368 !GOSUB 0 

360 !&X = LMC 

3 10 ! GOSUB 0 

371 !&X=ST 

3/2! GOSUB O 

37 3 ! &X=LST 

374! GOSUB 0 

375!&X=LO 

376JGOSUB 0 

37/ ! &X=LLO 

378 !GOSUB 0 

3 /Q ! &X=CF() 

380 !G()SUB 0 

381 ! &X=LCFO 
382 ! GOSUB 0 
383!&X=FL 
384! GOSUB 0 
385 !&X=LFL 
386 !GOSUB 0 
387 !&X=PS 
388! GOSUB 0 
38Q ! &X=LPS 
300! GOSUB 0 
30 I ! AX=02 
30 2! GOSUB 0 
303!AX=L02 
304 ! GOSUB 0 
3<35!&X=TF 
306 !GOSUB 0 
30 / ! & X=LTF 
3081GOSUB 0 
300! &X=TA 
400! GOSUB 0 
401 ! &X=LTA 
402 '.GOSUB O 
403 ! &X=MF 
404 ! GOSUB 0 
405!&X=LMF 
406! GOSUB 0 
407 ! &X=MA 
408 ! GOSUB 0 
400 ! & X=LM A 

4 10! GOSUB 0 

4 I I ! AX=WRT 
412! GOSUB 0 
4 l3!iX=LWRT 
4 1 4 ! GOSUB 0 
415! &X=EL 

4 l 6 ! GOSUB O 



4 | 7 ! &X=LEL 

4 1 8! GOSUB 0 

4 I 9 1 &X=FW 

420! GOSUB 0 

421 !&X=LFW 

422 IGOSUB 0 

423 ! &X=HLO 

424 IGOSUB 0 

425 ! &X=LRLO 

426 '.GOSUB 0 

427! &X=FBTF 

428 IGOSUB 0 

42P!&X=LFBTF 

430!GOSUB O 

431 !&X=FBTA 

432!GOSUB 0 

433 !&X=LFBTA 

434 IGOSUB 0 

435 ! &X=VFO 

436! GOSUB 0 

437 ! &X=LVFO 

4 38 IGOSUB 0 

439!&X=CFT 

440 IGOSUB 0 

44| ! &X=LFT 

442!G()SUB 0 

443 ! &X=CAUX 

444!G()SUB 0 

445 ! &X=LAU X 

446 '.G0SU3 0 

447 !&X=CAT 

448 IGOSUB 0 

44P !&X=LAT 

4 50 '.GOSUB 0 

451 ! PLOT DOT SCALE .5 

452 ! EX IT PART F 

453 !G()TO END 

454! #0 

4551WRITEF 3. AX 
456 1RTNSUB 
45 7 ! #END 

4581RUN NEW MIT. POLYFILE 

45 R ! END 
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MIT.ABCD. POLYFILE 
5- 7-84 I9i22:03 FUTIL 6.21 



OPENR 1 , "MIT. PASS. POLYDAT" 

OPENW 2 • " M IT. OUT . P()LY" 

PRNT OUTPUT FOR THE POLYGON MODULE WILL BE IN 'M IT . OUT . POLY' , 

at= "Polygon summary output file" 

WRITEF 2 ♦ AT 
AP=" 

WRITEF 2.AP 
WRITEF 2 , AP 

AT=" ITEM LCG WEIGHT" 



lOtWRITEF 2. AT 
1 1 'WRITEF 2 ,AP 
1 2 ! RE ADF 1 ,AX 
13! READF 1 , A Y 
1 4 ! &T- i 'CREW AND EFFECTS 


J '+AY + J ' 


J, + &X 


1 5 ! WRITEF 2. AT 
1 6 ! GOSUB RD 

1 7!AT="SLBM'S OR COMP. WATER 


"+AY+" 


»+&x 


1 8 ! WRITEF 2, AT 

1 Q ! GOSUB RD 

20 ! AT="S AN I TA RY TANKS 


" +A Y + " 


" + &X 


21 IWRITEF 2, AT 
22! GOSUB RD 

23 ! AT =i, LUBE OIL IN SUMPS 


"♦AY*" 


J, + &X 


24 ! WR ITEF 2. AT 
25 'GOSUB RD 

26 ! AT= J, F I XED CLEAN FUEL OIL 


"+AY+" 


" + &X 


2 7 ! W RITEF 2, AT 
28 'GOSUB RD 

29* AT®"**** TOTAL FIXED ITEMS * *** 


"♦AY+" 


"+&X 


30 ! WR ITEF 2 , AT 
3 1 ! GOSUB RD 

32! AT = "PROV ISIONS AND STORES 


J, +AY + J ' 


»+&x 


33 ! WR ITEF 2. AT 
34 'GOSUB RD 

35! AT="REVITALIZATlON 02 


"+AY+" 


" + &X 


36 ! WR ITEF 2, AT 
3 7 ! GOSUB RD 

38 ! AT="TORPEDOS IN FWD ROOM 


"+AY+" 


"+&X 


39 ! WRITEF 2 * AT 
40 ! GOSUB RD 

4 1 !AT="T()RPEDoS IN AFT ROOM 


"+AY + " 


"+&X 


42 ! WR ITEF 2, AT 
43! GOSUB RD 

44 ! AT = "TACT I CAL MISSILES FWD 


"+AY+" 


" + KX 


45 ! WR ITEF 2, AT 






40! GOSUB RD 

4 / ! &T=" TACT ICAL MISSILES AFT 


'• + A Y+" 


"+A.X 


48 ! WR I TEF 2, AT 
49! GOSUB RD 
50 ) AT =" WRT TANKS 


"♦AY* 41 


"+&X 


51 ! WRITEF 2, AT 
52! GOSUB RD 

o3!AT="4ESERVE ELECTROLYTE 


"+AY+" 


»+<<x 


54 ! WR ITEF 2 .AT 
55! GOSUB RD 

56 ! AT=" TOTAL FRESH WATER 


"+AY+" 


"+A X 
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5 / 1 WR ITEF 2. AT 
58 ! GOSlJB RD 








59! AT=“RESERVE LUBE OIL 
60 I WR ITEF 2. AT 
61 IGOSUB RD 


"+AY+ M 


M + A X 




62 ! AT="FUEL BALLAST FWD 
63 ! WR ITEF 2 , AT 
64 I GOSUB RD 


"+AY+" 


"+AX 




65 ! AT = i 'FU EL BALLAST AFT 
66 ! WR ITEF 2 .AT 
67 ! GOSUB RD 


<, **'AY+ 4 ' 


<l + AX 




od ! AT= M VAR I ABLE FUEL OIL 
69 I WR ITEF 2, AT 
70 IGOSUB RD 
/ 1 I WR ITEF 2 , AP 
72 I WR ITEF 2 • AP 

/3! AT= M VARI ABLE BALLAST TANK 
74 ! WR ITEF 2 • AT 


"♦AY*” 
DATA * M 


"+AX 




75 I AT=" TANK GROUP 

76! WRITEF 2, AT 


LCG 


CAPACITY (TONS)" 


77 ! AT ="FORW ARD TRIM 
7d I WRITEF 2, AT 
79 IGOSUB RD 


" +A Y+ " 




11 +AX 


dO I AT= M AUX IL I ARY 
81 I WR ITEF 2, AT 
82 I GOSUB RD 


» +A Y+ " 




" +AX 


83 I AT="AFTER TRIM 

84 I WR ITEF 2, AT 

85 i G oto end 
86i#RD CONTINUE 
87 I READF 1 ,AX 
88IREADF I , AY 
891 RTNSUB 
90!*END 
91 I END 


"+AY+’*' 




J, +AX 



Thesis 
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Thesis 
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